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20.  ABSTRACT  (continued) 
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to  monitor  the  presence  of  excess  oxygen  in  silicon  by  the  presence 
of  this  peak  in  the  recombination  spectra  of  silicon  after  irradiation. 

Thus,  observation  of  0.790  eV  luminescence  from  float  zone  grown 
silicon  samples  implanted  with  arsenic  through  silicon  dioxide,  and 
o»;t  from  identical  samples  in  which  the  silicon  dioxide  was  chemically 
removed  prior  to  arsenic  implantation,  is  interpreted  as  direct  evidence 
* ' for  the  recoil  implantation  of  oxygen  into  the  silicon  from  the  silicon 

dioxide  layer.  The  effect  of  recoil-implanted  oxygen  on  the  transport 
propertj.es  of  silicon  implanted  with  arsenic  to  fluences  of  10^5(-ni"2 
and'5  x lO^^cm"^  is  then  examined.  Sheet  resistivity  and  sheet  Hall 
coefficient  measurements  following  various  annealing  stages  indicate 
identical  recovery  for  through-oxide  implants  and  for  bare  surface 
implants.  ,^Comparison  of  the  impurity  profiles  obtained  by  secondary 
ion  mass  spectroscopy  with  free  carrier  profiles  obtained  by  differential 
Hall  effect  measurements  shows  identical  activation  for  through  oxide 
and  bare  surface  arsenic  implants  for  both  600°C  and  1000°C  anneals. 
Additionally,  mobility  profiles  obtained  during  the  differential  Hall 
effect  measurements  indicate  mobility  values  for  the  bare  surface  and 
through  oxide  implants  that  agree  with  values  expected  for  arsenic  doped 
silicon.  Measurements  obtained  in  this  study  are  consistent  with  the 
interpretation  that  residual  damage  seen  for  through  oxide  implanted  samples 
is  caused  by  random  nucleation  during  the  recrystallization  of  the  implanta- 
tion-induced amorphous  layer,  due  to  a reduction  in  the  epitaxial  regrowth 
rate  of  the  implanted  layer  by  the  presence  of  excess  recoil-implanted 
oxygen . 
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The  properties  of  silicon  after  annealing  of  arsenic  implantations 
through  silicon  dioxide  are  studied,  and  compared  with  those  resulting  from 
arsenic  implantations  into  bare  silicon  surfaces.  Photoluminescence  studies 
of  radiation  damaged  silicon  establish  the  0.970  eV  peak  as  arising  from  a 
carbon-dependent  center.  The  same  studies  also  confirm  that  the  0.790  eV 
luminescence  arises  from  an  oxygen  dependent  center,  and  suggest  that  this 
center  may  depend  on  carbon  as  well.  With  the  oxygen  dependence  of  the  center 
responsible  for  the  0.790  eV  luminescence  established,  it  is  possible  to  monitor 
the  presence  of  excess  oxygen  in  silicon  by  the  presence  of  this  peak  in  the 
recombination  spectra  of  silicon  after  irradiation.  Thus,  observation  of 
0.790  eV  luminescence  from  float  zone  grown  silicon  samples  implanted  with 
arsenic  through  silicon  dioxide,  and  not  from  identical  samples  in  which  the 
silicon  dioxide  was  chemically  removed  prior  to  arsenic  implantation,  is  inter- 
preted as  direct  evidence  for  the  recoil  implantation  of  oxygen  into  the 
silicon  from  the  silicon  dioxide  layer.  The  effect  of  recoil-implanted  oxygen 

on  the  transport  properties  of  silicon  implanted  with  arsenic  to  fluences  of 
15  -2  15  -2 

10  cm  and  5 x 10  cm  is  then  examined.  Sheet  resistivity  and  sheet  Hail 
coefficient  measurements  following  various  annealing  stages  indicate  identical 
recovery  for  through-oxide  implants  and  for  bare  surface  implants.  Comparison 
of  the  impurity  profiles  obtained  by  secondary  ion  mass  spectroscopy  with  free 
carrier  profiles  obtained  by  differential  Hall  effect  measurements  shows 
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identical  activation  for  through  oxide  and  bare  surface  arsenic  implants  for 
both  600°C  and  1000°C  anneals.  Additionally,  mobility  profiles  obtained 
during  the  differential  Hall  effect  measurements  Indicate  mobility  values 
for  the  bare  surface  and  through  oxide  implants  that  agree  with  values 
expected  for  arsenic  doped  silicon.  The  measurements  obtained  in  this  study 
are  consistent  with  the  interpretation  that  residual  damage  seen  for  through 
oxide  implanted  samples  is  caused  by  random  nucleation  during  the  recrystal- 
lization of  the  implantation-induced  amorphous  layer,  due  to  a reduction  in 
the  epitaxial  regrowth  rate  of  the  implanted  layer  by  the  presence  of  excess 
recoil-implanted  oxygen. 
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1.  INTRODUCTION 

Ion  implantation  has  evolved  from  an  area  of  specialized  research 

investigations  into  a technique  routinely  used  throughout  the  semiconductor 

industry.  Ion  implantation  as  a process  tool  has  many  advantages  over 

diffusion  technology,  including  more  accurate  doping  control  and  better 

lateral  registration,  lower  temperature  processing,  better  control  of  the 

depth  and  profile  of  the  dopant  impurities,  and  the  potential  for  novel 

d- .icc  structures  [1,2],  The  implantation  dose  can  be  accurately  controlled 
16  "2 

from  about  10  cm  , which  is  equivalent  in  impurity  concentration  to  the 

11  -2 

heaviest  diffused  predeposition,  down  to  a dose  of  10  cm  , which  has  no 
counterpart  in  diffusion  technology  [2]. 

Silicon  device  technology  provides  an  excellent  illustration  of  the 
advantages  of  ion  implantation.  For  many  applications,  arsenic  would  be  a 
more  desirable  dopant  impurity  than  phosphorus,  due  in  large  part  to  its 
lower  diffusion  coefficient  [3],  and  hence  lesser  redistribution  during 
subsequent  processing  steps.  Unfortunately,  arsenic  diffusion  technology  is 
difficult  [2],  while  arsenic  implantation  is  fairly  simple.  High-dose 
arsenic  implantations  are  commonly  used  to  form  buried  layers,  emitters  for 
bipolar  transistors,  and  for  doping  of  contact  regions  [1,2],  making  arsenic 
the  donor  impurity  of  choice  for  many  applications  which  were  unrealistic 
prior  to  the  advent  of  ion  implantation. 

Despite  its  many  advantages,  ion  implantation  would  not  have 


achieved  its  widespread  popularity  were  it  not  for  its  excellent  com.patibi lity 
with  the  silicon  planar  process.  The  use  of  silicon  dioxide  layers  for  surface 
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passivation  and  selective  area  masking  is  an  essential  part  of  the  silicon 
planar  process  [IJ.  Therefore,  the  discovery  that  arsenic  implantations 
through  the  edges  of  tapered  oxide  layers  lead  to  densely  packed  defect 
structures  which  cannot  be  removed  by  high  temperature  annealing  [4J  is  a 
matter  of  great  concern.  In  the  experiment  that  led  to  this  discovery, 
fluences  of  lo''^  to  10^^As"^/cm^  were  implanted  into  silicon  substrates 
partially  masked  by  islands  of  tapered  oxide,  and  then  the  samples  were 
annealed  for  30  minutes  at  1000°  C to  simulate  the  growth  of  an  epitaxial 
layer.  Examination  of  the  samples  by  t'ransmission  electron  microscopy  then 
revealed  a stable,  high  defect  density  region  at  the  tapered  edges  of  the 
oxide  mask.  Significantly,  the  annealing  temperature  used  is  considerably 
higher  than  the  600°  C annealing  normally  required  to  restore  the  crystal- 
linity of  heavily  implanted  layers  [1,.  In  the  same  study  it  was  also 
discovered  that  identical  disorder  could  be  reproduced  by  implantations  through 
uniformly  thin  oxides.  In  both  cases,  no  residual  damage  could  be  seen  in 
regions  from  which  the  silicon  dioxide  was  removed  prior  to  implantation. 
Because  identical  damage  was  seen  for  uniformly  thin  cxide  layers,  stress 
effects  at  the  interface  bettireen  the  bare  surface  and  the  cxide  edge  could 
be  ruled  out.  Instead,  the  damage  was  produced  by  the  penetration  cf  the 
arsenic  beam  through  the  oxide  hayer.  As  the  arsenic  icxis  penetrate 
the  silicon  dioxide  layer,  they  suffer  colli.sions  with  oxygen  atoms  and  thus 
could  potentially  knock  large  concentrations  of  oxygen  into  the  silicon 
substrate.  Tliis  aspect  of  the  problem  was  examined  by  Chu  and  co-workers  [3] 
using  Che  c-’chnique  of  helium  ion  backscattering  to  measure  the  extent  of 
residual  disorder.  Tliey  discovered  that  implantations  of  arsenic  and  oxygen 
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into  bare  silicon  surfaces  reproduced  the  damage  seen  from  arsenic  implan- 
tations through  thin  oxide  layers,  and  that  arsenic  implantations  into  bare 
silicon  surfaces  also  lead  to  residual  disorder  [6].  The  problem  was  later 
examined  by  transmission  electron  microscopy  [7],  which  showed  that  arsenic 
implantations  into  bare  silicon  produced  general  lattice  disorder,  but 
that  this  disorder  was  significantly  less  and  defects  were  more  widely 
spaced  than  when  both  arsenic  and  oxygen  were  implanted.  This  study  also 
showed  that  arsenic  and  oxygen  implantations  into  bare  silicon  yielded  the 
same  type  of  damage  as  seen  from  arsenic  implantations  through  silicon 
dioxide.  Thus,  while  these  results  implicate  recoiling  oxygen  atoms  as 
being  responsible  for  the  residual  damage,  it  has  not  been  possible  to 
identify  the  presence  of  oxygen  recoils  in  the  through-oxide  implanted  layers. 
While  arsenic  has  been  found  to  interact  strongly  with  oxygen  impurities  in 
the  annealing  of  radiation-damaged  silicon  [8] , none  of  these  studies  have 
examined  the  role  of  arsenic  in  the  formation  of  the  residual  defects.  Finally, 
and  of  greatest  importance,  none  of  these  studies  have  examined  the  effects 
of  these  residual  defects  on  the  electrical  properties  of  the  implanted  layers. 
The  aim  of  the  present  work  is  to  examine  the  problems  created  by  through- 
oxide  arsenic  implantation,  and  to  contrast  the  properties  of  these  layers 
with  those  of  arsenic  implantations  into  bare  silicon  surfaces,  to  determine 
the  causes  of  this  residual  damage  and  the  extent  of  its  effects. 
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2.  RECOIL  IMPLANTATION 

Much  of  the  work  reported  here  involves  the  interactions  that  occur 
when  an  energetic  ion  strikes  a solid  target.  Analysis  of  these  events  leader 
to  predictions  of  the  depth  an  energetic  ion  will  penetrate  into  the  solid, 
and  the  extent  of  damage  created  by  transfer  of  the  ion's  incident  energy 
to  the  target  atoms.  It  is  also  possible  to  estimate  the  energy  and  flux 
of  recoiling  target  atoms,  and  to  predict  the  erosion  of  target  surfaces  due 
to  sputtering.  This  chapter  provides  a review  of  the  assumptions  and  approxi- 
mations needed  to  describe  these  processes,  with  specific  appi;i-.u  lor  ^o  t.!.e 
case  of  As  implantation  through  Si02  layers  into  silicon  substrates.  The 
comprehensive  treatment  by  Sigmund  [9]  of  the  Lindhard,  Scharff,  Schi^tt  (LSS) 
theory  [10]  forms  the  basis  for  this  description.  The  theory  is  presently 
incomplete  and  is  somewhat  controversial,  especially  in  regard  to  implantations 
into  multilayer  targets.  Therefore  our  aim  here  will  be  to  describe  in  general 
terms  the  physical  mechanisms  that  occur  and  to  develop  reasonable  estimates 
of  ion  and  recoil  behavior. 

As  an  incident  ion  strikes  a solid  target  and  eventually  comes  to 
rest,  it  deposits  its  kinetic  energy  into  the  target  through  two  mechanisms: 

1)  transfer  of  kinetic  energy  to  target  atoms  in  essentially  elastic  collisions 
(referred  to  as  nuclear  stopping);  and  2)  transfer  of  energy  in  the  excitation 
of  the  orbital  electrons  of  the  target  atoms,  ions,  or  both  (electronic  stopping). 
These  mechanisms  occur  simultaneously  but  are  treated  as  being  distinct  for 
ease  in  analysis.  Nuclear  stopping,  being  conceptually  the  simpler  of  the 
two,  will  be  treated  first.  We  consider  the  ion-target  collision  as  the 
classical  problem  of  the  ion's  motion  in  the  central  field  of  the  ion-atom 


potential  [9j. 


For  the  case  of  the  target  atom  initially  at  rest,  the  conservation 
of  energy  and  momentum  require  [11,12J 


tan  = 


M,  sin9 
Mf  +M2  sin9 


6"  = 


TT  - 9 


E'  = E - T 
E"  = T 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


2 9 

X = T sin  f 
m 2 

4M^M2 

^m  ~ ^ 


(2.5) 

(2.6) 


where  is  the  mass  of  r'le  incident  ion  and  is  the  target  atom  mass;  E is  the 
ion's  initial  encTf,/;  I is  Che  energy  transferred  in  the  collision;  ''  is 
the  scattered  angle  of  the  projectile  in  the  laboratory  system  with  energy 
E'  after  the  collision;  and  is  the  angle  taken  by  the  scattered  particle 
in  the  laboratory  system  with  energy  E"  after  the  collision,  wi-..  the 
scattering  angle  in  the  center  of  mass  system  (Fig.  2.1).  Here  electronic 
excitation  is  considered  only  a.s  a source  of  energy  loss,  and  is  assumed  to 
have  negligible  effect  on  the  collision  dynamics.  The  experience  from  ion- 
gas  collisions  in  the  keV  region  [13]  is  that  electronic  energy  loss  amounts  to 
only  a few  percent  of  the  energy  transfer  and  thus  is  assumed  to  have  a negli- 
gible effect  on  either  the  amount  of  energy  transfer  or  the  post -collision 
trajectories  of  the  ion  and  recoiling  target  atom. 
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Using  trigonometric  identities,  the  Eqs.  (2.2)  and  (2.5)  yield  a 


more  useful  form 


cos  i>" 


= (f)' 


giving  the  trajectory  of  the  recoil  in  the  laboratory  frame.  To  find  the 
angle  0,  we  integrate  the  potential  between  ion  and  target  atom  V(r)  | 11,12] 
in  the  center  of  mass  system 


9 = n - 2p 


R . '/i-v(r)|e  -p^(r„^ 

min  ' r 2 


where  p is  the  impact  parameter  and 

M^E 

^r  " (Mj^  i-M^) 

is  the  particle  energy  in  the  center  of  mass  system.  R . is  the  minimum 
distance  of  approach  (Fig.  2.1)  determined  by  the  zero  of  the  ra.li.ccnd  ;n 
the  denominator  of  the  integral  in  Eq . (2.8).  Inversion  of  Eq.  (2.8)  yields 
the  differential  nuclear  cross-section  do(9)  assuming  azimuthal  symmetry: 


da  = -2n  p ^ d9  = ^ dfi 

<39  sin  9 Q0 


(2.10) 


where 


dCl  = 2tt  sin  9 d9 


(2.11) 


Large  angle  collisions  are  a rarity  at  energies  greater  than  1 keV ; 
thus,  typical  impact  parameters  are  large.  For  these  cases  (the  majority 


of  all  stopping  events),  we  have  [9]  V(R)/E^  « 1 and  thus  expansion  In  terms 
of  V/E^  yields  [14] 


It  is  important  to  note  that  this  new  expression  is  linear  in  V(R);  thus, 
simultaneous  collisions  can  be  correctly  described  by  the  theory  in  this 
approxima  tion . 

To  perform  useful  calculations,  one  must  find  an  expression  for 
the  interatomic  potential  V(r).  Since  accurate  empirical  potentials  are 
available  for  only  a few  ion-target  combinations,  theoretical  models  are 
employed.  The  Thomas-Fermi  model  is  most  often  used,  due  to  its  simplicity 
and  generality.  This  model  essentially  considers  the  ion-target  system  as 
two  Coulombic  nuclei  surrounded  by  an  electron  gas  of  variable  density, 
considered  uniform  on  a microscopic  scale.  From  this  picture  of  the  atomic 
system,  minimization  of  the  total  energy  by  a variational  procedure  leads  to 
the  Thomas-Fermi  nonlinear  differential  equation  [9,15). 


V^(V-V  ) = Anu  e(V-V  (2.13) 

^ o b ' o ' 

where  V is  the  electrostatic  potential  of  the  system,  'cV^  is  '■he  Fenri 

energy  , and  is  given  by 

a = (2.1A) 

o 

where  a is  the  screening  radius  chosen  by  Lindhard  ^ al^.  [lOJ 

2/3  2/3  -1/2 

a = 0.8853  a (Z , ' -I- Z,  ' ) (2.15) 

o i z 

with  Zj^  and  the  atomic  numbers  of  the  ion  and  target  nuclei,  respectively, 
and  a^  the  Bohr  radius  (0.528  X). 


I 

[ 
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The  solution  to  this  equation  is  then 

V(r)  = fg(R/a) 


(2.16) 


where  4^  is  a function  found  by  numerical  integration.  Conveniently,  the 
function  4^  can  be  approximated  by  inverse  powers  of  its  argument  over 
extended  distances  (Fig.  2.2).  The  form  of  Eq . (2.16)  illustrates  the 
utility  of  the  Thomas-Fermi  approach,  in  that  the  details  of  atomic 
structure  do  not  enter  into  the  calculations,  and  thus  the  theory  can  easily 
be  extended  to  cover  all  possible  ion-target  combinations.  The  only 
changes  that  occur  from  one  ion  target  system  to  another  are  in  length  scale 
(given  by  Eq . (2.15)),  and  the  total  number  of  electrons  (determined  by 
and  Z^) . Using  Eq . (2.16)  for  the  interatomic  potential,  Eq . (2.8) 


becomes 


£ f 


d(R/a)/(R/a) 


9 = -n  - 2 - j 

3 , 

Rmrn 


^1^2  ^(R/a)  (P/®) 


(2.17) 


1 - 


aE^  R/a 


(R/a)' 


iuLS  expression  tan  be  simplified  by  a;. 
ly  lo  erst  power;  c'i  Its  atgurtnr  (F  u.  2.2) 


4^  (R/a)  = --  (a/R)®*’-  (2.18) 

with  s an  integer  and  a numerical  constant.  Ihe  expression  for  9 then 
results  by  substituting  Eq . (2.18)  into  the  approximate  form  Eq . (2.12) 

V k , 

® ® (2. 19) 


r 

f 


f 
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LP-1270 


Figure  2.2.  Thomas-Fermi  screening  function  5q(x)  for  neutral 

atoms  and  its  approximation  by  inverse  powers  of  its 
argument  (after  Gombas  [16]). 
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where 


1 3s-l 

Y =»  - 

S 2 

and  € is  a dimensionless  scaling  unit  for  the  center  of  mass  energy 

aE 

€ . £— 


Z^Z26 


(2.20) 


(2.21) 


Having  finally  obtained  an  expression  for  9,  we  now  can  substitute  into 
Eq . (2.10)  to  obtain 

dCT(0)  = (constant)  9’^  d9  (2.22) 

A more  useful  expression  is  obtained  by  substituting  for  9 from  Eqs.  (2.5)  and 


(2.6)  assuming  9 is  small: 


da(T)  = dT 


(2.23) 


with 


and 


1 

m = — 
s 


(2.24) 


m „ 2 

2 (^1^  /^^1^2®  ' 


C = ^ \ { ) 2m 

2 m \M„/  V a / 


(2.25) 


Fitting  constants  to  the  Thomas-Fermi  function  yields  [16] 


S/3  = '-30^ 


S/2  = 0-327 


= 0.5 


(2.26) 

(2.27) 

(2.28) 


J 


r 


"1 


The  approximations  used  for  various  energy  intervals  are: 
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m = 

1/3 

for 

6 < 0.2 

(2.29) 

m = 

1/2 

for 

0.08  ^ t <0.2 

(2.30) 

m - 

1 

for 

€ > 10 

(2.31) 

To  this  point  we  have  considered  a as  a function  of  two  variables, 

€ and  9 (or  equivalently  E and  T).  However,  Lindhard  e^  a_l.  have  shown  [10] 
that  to  the  validity  of  the  power  approximation,  Eq . (2.18),  the  1 f fcrci'.tia  1 
cross  section  can  be  considered  a function  of  the  single  variable  t,  where 


t = 6 


2 T 


(2.32) 


Thus  Lindhard  e^  aj^.  [10]  use 


da  = na^  f(t^^^) 


2t 


3/2 


(2.33) 


1/2, 


where  f(t  ) is  plotted  in  Fig.  2.3.  The  next  step  is  to  generalize  from 
a collision  between  an  ion  and  a single  target  atom  to  the  case  of  a target 
consisting  of  a collection  of  atoms  distributed  randomly  with  a density  N 

3 

(atoms/cm  ) over  a depth  Ax.  Fcr  this  case,  the  probability  P for  an  ion  to 
suffer  a collision  specified  by  a cross  section  do  is  given  by 

dp  = nAx  do  (2.34) 

for  Ax  small  enough  so 'that  dP  « 1.  For  an  initial  energy  E with  recoil  energy 
T',  tile  average  energy  loss  VAK)  is  then 


(AE>  = Tl'dP  = NAx  ''  T ^ dT 
■ T=0  ' ^ 


(2.35) 


In  the  limit  of  infinitesimal  Ax,  Eq . (2.35)  becomes 


Figure  2.3.  Differential  cross  section  for  elastic  collisions 
^ 1 _2 
f(t‘^)  vs.  t‘  ; with  approximation  by  r potential 

and  high  energy  limit  of  Rutherford  scattering 

(after  LSS  [10]). 
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Figure  2.4.  Nuclear  stopping  (S^(€)),  electronic  stopping 

(Sg(€))>  and  total  stopping  in  dimension- 

less units.  The  electronic  stopping  power  is  given 

for  k = 0.12.  Also  shown  is  the  approximation  of 

_2 

the  nuclear  stopping  power  for  a r potential  s °(€) 

n 

(After  LSS  [10] . ) 
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Thus  we  have  finally  arrived  at  an  analytic  description  of  the 
elastic  energy  transfer  during  the  ion-target  collision.  Naturally,  this  is 
the  energy  used  to  produce  recoiling  target  atoms,  and  thus  will  be  most 
important  for  estimation  of  recoil  densities.  The  development  of  these 
expressions  is  rigorous  within  the  assumptions  made,  and  is  firmly  based  on 
classical  dynamics.  It  is  important  to  note  that  quantum  effects  (i.e., 
deviations  from  classical  mechanics)  do  indeed  .take  place.  These  effects 
include  purely  electronic  excitation,  exchange  scattering,  and  resonant 
charge  transfer  [17],  however,  these  effects  typically  amount  to  only  a few 
percent  of  the  energy  exchange  for  ions  with  kinetic  energies  of  1 keV  or 
greater. 

We  now  consider  electronic  stopping  for  the  case  of  implanted  ions. 
For  these  conditions,  the  majority  of  target  electrons  are  travelling  much 
faster  than  the  incident  ion,  so  that  the  electrons  cannot  pick  up  any  energy 
from  direct  collisions  with  the  ton. 

Earlier  in  this  discussion  it  was  pointed  out  that  electronic 
excitation  amounts  to  only  a few  percent  of  the  energy  transfer  during  nuclear 
collisions.  However,  the  major  contribution  to  electronic  stopping  occurs 
from  chose  collisior.s  in  which  the  nuclei  are  far  enough  apart  to  prevent  any 
appreciable  transfer  of  kinetic  energy,  but  in  which  some  overlap  of  outer 
shell  electrons  can  occur  [9].  Since  the  differential  cross  section  is  propor 
tional  to  the  impact  parameter  (Eq . (2.10)),  these  collisions  can  assume  a 
considerable  statistical  weight.  There  is  no  definitive  theory  on  electronic 
stopping;  therefore,  we  again  concentrate  on  understanding  die  physical 
mechanisms  rather  dian  focusing  on  details  of  the  computational  procedure. 
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In  the  coordinate  system  of  the  ion,  orbital  electrons  have  a small 
drift  velocity.  llius , when  an  electron  is  scattered  by  the  moving  ion  [10], 
or  is  captured  by  the  moving  ion  [18],  a net  momentum  (which  is  proportional 
to  the  ion  velocity)  is  transferred  to  this  electron.  According  to  these 
arguments,  the  energy  lost  by  the  ion  would  then  yield  a stopping  power 
proportional  to  the  ion's  velocity 


S = kv  = kE 
e 


(2.42) 


The  value  for  the  constant  k depends  on  the  model  used;  however,  the  common 
practice  in  compilations  of  range  statistics  is  to  use  experimentally- 
determined  stopping  coefficients  when  available  [19]. 

It  should  also  be  noted  that  when  Hartree-Foch  models  are  used  for 
the  ion-target  combination,  pronounced  oscillations  are  found  for  the 
electronic  stopping  power  as  a function  of  the  ion  atomic  number  [20|.  Such 
oscillations  are  not  predicted  by  the  simpler  models  mentioned  above,  but 
are  observed  exper im.enta lly  [21j  and  have  been  included  in  compilations  of 
range  statistics  [19].  Recent  w’ork  [22]  has  shower.,  however,  that  their  effects 
have  been  overestimated  for  random  implants  into  silicon  targets. 

Having  developed  expressions  for  the  nuclear  and  electronic  stopping 
powers,  we  are  now  in  a position  to  formulate  the  depth  distribution  of 
implanted  ions.  Along  an  incremental  path  length  Ax,  an  incident  ion  suffers 
an  energy  loss  due  to  nuclear  and  electronic  stopping,  given  by  [10]: 


Ae  = - N(S  (E)  + S (E))Ax 
n e 


(2.43) 
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In  the  limit  of  vanishingly  small  Ax 


^ - N(S  (K)  + S (E)) 

dx  n e 


(2.44) 


so  that  to  a first  approximation 


0 j € 

" ,_dlL_  = '>  _dE 

- dE/dx  “ N(S^(E) 


(2.45) 


or  in  €-0  units 


”<e>  • ^ sjo 


(2.46) 


Unfortunately,  the  range  travelled  by  an  ion  is  not  a physical  observable, 
since  the  ion  does  not  travel  in  a straight  line  from  the  target  surface. 
Instead,  the  ion  is  deflected  at  an  angle  f each  time  it  suffers  a nuclear 
collision  (Fig.  2.1a).  Because  the  ions  approach  target  atoms  with  a variety 
of  impact  parameters,  we  obtain  a distribution  of  ion  penetration  depths, 
or  equivalently,  a probability  distribution  for  the  depths  a particular  ion 
can  come  to  rest  within  the  sample  -surface. 

There  are  several  approaches  to  the  calculations  of  ion  ranges, 
but  none  of  these  methods  produce  less  than  a 10  - 20%  deviation  from  measured 
values  [9|.  Here  we  follow  the  method  of  Sigmund  [9],  which  most  easily  allows 
the  calculation  of  recoil  densities.  Tlie  method  of  Brice  has  also  recently 
been  modified  to  include  the  effects  of  recoiling  target  atoms  [23].  Sigmund 
employs  the  concept  of  a final  range  distribution  F(E,e,r),  which  is  the 
probability  that  an  ion  with  initial  energy  E directed  along  a unit  vector  e 
will  come  to  rest  a (vector)  distance  from  the  point  at  which  it  encountered 
the  sample  surface  (r  = 0).  Fairly  straightforward  arguments  [9]  yield  the 


equation 
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^ — f*  , — 

- cost  ^ F(E,e,r)  = NJdal  F(E  ,e  ,x)  - F(E-T,e'x)] 


+ i^F^(E,c-,x) 

where  x is  the  depth  measured  normal  to  the  target  surface  and 

cos  t = e -x  ( 

with  e'  a unit  vector  in  the  direction  of  the  scattered  ion,  and  x a unit 
vector  in  the  x direction.  Here  we  have  already  integrated  the  y and  Z 
dependence.  This  equation  is  basically  a Taylor's  series  expansion, 
valid  for  terms  linear  in  5x.  For  symmetry  reasons  [9],  we  have  for 
arbitrary  incident  direction  : 

F(E,e,x^  = F(E,cos  'y , x)  ( 


together  with  the  normalization 


/dx  F (E  , cos  t > x)  = 1 


where  the  integration  extends  over  all  space  to  include  backscattered  ions. 
We  now  have  a compatible  set  of  equations  which  can  ideally  be  solved  for 
the  di.stribution  of  penetration  depths. 

Most  calculations  do  not  attempt  to  solve  Eqs.  (2.47)  and  (2.50). 
but  rather  expand  these  equations  in  terms  of  spatial  moment  F^(E,e)  over 
the  depth  distribution  F(E,e,r)  [9,24]; 


F*^(E,e)  = Jdx  x'^  F(E  ,e  ,x) 


with  F =1.  This  procedure  yields 
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nt'^'^E.e)  = Nj'da{K^(E,e)  ' F’^(E-T,e)} 

+ N S^(E)  ^ F"(E,e)  (2.52) 

e sE 

wi th  F°  = 1 . 

While  the  expansion,  Eq . (2.46),  relieves  many  computational 

dif liculties , it  must  be  noted  that  Eq . (2.46)  requires  knowledge  of  the 

previous  F^(E)  to  determine  F^'"^^(E).  Tnis  places  great  demands  on  the 

accuracy  to  which  the  previous  coefficient  mus*:  be  obtained.  It  can  be 

shown  that  the  first  moment  yields  the  average  penetration  depth  [9]  and 

a Gaussian-type  distribution  centered  around  R with  a standard  deviation 

P 

or  straggle  (referred  to  as  ^p)  • many  practical  applications,  the 

gaussian  apprcximation  is  sufficient.  For  example,  Marcyk  and  Streetman  [25] 
have  shown  that  Che  gaussian  approximation  describes  the  summation  of 
implantations  at  various  energies  into  a common  Si  target. 

Thus  far  v;e  have  described  range  distributions  into  an  amorphous 
substrate  of  only  one  chemlca''  species.  To  include  cuompound  substrates,  a 
weighted  average  tf  :he  stopping  powers  of  each  elemetit  over  their  elemental 
density  is  performed.  Thus  Si02  is  modelled  by  a p-seudo-element  of  atomic 
Tiamher  10  and  atomic  weight  20.0.  An  entirely  different  effect  occurs  for 
ions  incident  cn  single  crystalline  targets.  Tlii.s  effect,  referred  to  as 
''channelling"  [26],  occurs  when  an  icT.  is  incident  upon  (or  scattered  into)  a 
path  within  a few  degrees  of  a low  index  crystalline  direction.  For  these 
ions,  a collective  steering  effect  from  rows  of  target  atoms  repels  the 
channelled  ions  and  prevents  them  from  suffering  direct  nuclear  collisions. 
These  ions  thus  travel  further  due  to  their  more  direct  path,  encountering 


only  electronic  stopping.  Since  a randomly-oriented  silicon  substrate  is 
driven  amorphous  by  doses  of  ~ 10^^  ions/cm  ^ [27],  it  is  possible 
that  a significant  fraction  of  the  initial  dose  could  become  channelled 
before  the  implanted  silicon  undergoes  an  amorphous  transition  [28].  This 
phenomenon  would  then  yield  a tail  near  the  end  of  the  distribution.  Since 
many  of  the  properties  of  semiconductor  devices,  such  as  junction  capacitance 
[28]  and  solar  cell  collection  efficiences  [29],  are  sensitive  to  the  shape 
of  dopant  distributions,  channelling  may  be  an  important  consideration. 

Simple  methods  of  controlling  or  eliminating  this  effect  could  be  of  great 
technological  importance. 

We  finally  consider  the  problem  of  recoil  implantation.  As  the 
incident  ion  encounters  target  atoms,  it  transfers  some  of  its  kinetic  energy 
to  the  substrate  atoms  according  to  Eqs . (2.2),  (2.4),  and  (2.5).  At  somewhat 
lower  energies,  the  incident  ion  continues  in  a modified  direction.  Should 
the  energy  transfer  be  greater  than  the  chemical  energy  binding  the  target 
atoms  in  place,  the  target  atoms  would  proceed  in  the  direction  specified  by  i' 
(determined  by  Eq . (2.7)),  These  atoms  travel  a distance  which  depends  on 
their  energy  loss  in  the  target  medium.  Such  primary  recoils  might  then  set 
other  target  atoms  in  motion  due  to  further  nuclear  collisions,  and  sc  on. 

To  gain  an  understanding  of  the  distribution  of  recoil-implanted 
atoms,  one  first  considers  the  case  of  a densely-packed  monomolecular  layer  on 
the  surface  of  a silicon  target.  Such  a situation  could  be  approximated  by 
hydrocarbon  pump  oil  adsorbed  onto  a wafer  surface  in  the  target  chamber  of  an 
ion  imp’anter.  Since  the  layer  is  only  one  molecule  thick,  electronic  stopping 


can  be  neglected,  and  hence  we  only  consider  nuclear  collisions. 
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Considering  an  As  beam  of  100  keV  striking  a C layer,  we  see  from 
Eq.  (2.6)  that  at  most  48  keV  of  energy  could  be  transferred  to  a carbon 
atom,  for  a direct  collision  (p  = 0).  This  would  yield  a penetration  depth 
for  the  carbon  atom  of  1118  K,  as  estimated  from  projected  range  compilaticr.s 
[19J,  whereas  the  arsenic  ion  would  only  have  a projected  range  of  582  X. 
However,  this  condition  is  only  for  the  extremely  unlikely  event  of  a 
direct  collision;  more  probably  the  incident  beam  will  have  a variety  of 
impact  parameters,  which  would  yield  a distribution  of  recoil  energies  with 
accordingly  smaller  projected  ranges.  Additionally,  these  lower  energy 
recoils  will  travel  at  an  angle  to  the  surface,  determined  by  Eq.  (2.7),  and 
thus  would  have  shallower  penetration  depths  measured  normal  to  the  sample 
surface.  It  is  important  to  note  that  the  amount  of  energy  transfer  depends 
on  the  ion-atom,  mass  ratio.  For  the  case  of  a nitrogen  beam  striking  a carbon 
layer,  the  maximum  energy  transfer  could  be  as  large  as  99%. 

To  estimate  the  recoil  yield  from  a thin  target  into  a substrate  of 
a different  material,  Moline  e_t  al^.  [31]  have  developed  a simple  formalism 
to  estimate  recoil  implantation  yields  for  the  case  of  a heavy  projectile 
impinging  on  a lighter  target.  If  the  incident  ion  strikes  a target  atom  a 
distance  ? from  the  target-substrate  interface,  this  recoil  can  penetrate  the 
substrate  only  if  the  range  of  the  recoil  projected  normal  to  the  target 
surface  is  greater  than  the  distance  to  the  interface. 

Mathematically,  if  the  collision  transfers  an  energy  T to  the  target 
atom,  we  consider  the  recoil  only  for  the  case  where 


R (T)  cos  4"  > 5 
P 


(2.53) 
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The  total  number  of  recoils  n from  depth  is  then  given  by 


^ = daX 


(2.54) 


where  is  the  minimum  value  of  energy  transfer  for  which  Eq . (2.53)  holds, 
is  the  density  of  target  atoms,  is  defined  by  Eq . (2.6)  and  dC7(T)  is 
the  differential  cross  section  for  a collision  with  energy  transfer  T. 

Using  Eqs.  (2.23),  (2.25)  and  (2.30)  to  represent  da(T),  Eq . (2.54) 


becomes 


2 1/2 
, 0.33TTa  N 

dn t 1 /_m  I 

ds:  'Stj  ■ - 


(2.55) 


with  a defined  in  Eq . (2.15)  and  defined  in  Eq.  (2.21).  Ihe  use  of  Eq . (2.30) 

2 

implies  the  approximation  of  the  Thomas-Fermi  screening  function  by  a 1/r 
potential.  Moline  e^  al^.  further  approximate  the  stopping  powers  of  the 
target  for  both  the  incident  ion  (S^)  and  recoil  ($2)  by  constants  (Fig.  2.4). 
With  this  approximation,  we  have 


E^  = E - (W-=^)S^ 


(2.56) 


where  E is  the  incident  ion  energy,  E^  the  projectile  energy  at  distance 
from  the  target-substrate  interface,  and  W the  thickness  of  the  target 
layer  (i.e.  ' = W is  the  front  surface  of  the  target). 

Using  Eq . (2.56),  Eq . (2.54)  becomes 


n = 0.33"a^N^  I •t'  F' 
0 1 


(2.57) 


T 


where 


VE.,'" 

'57J  • ' '=1  » =2^ 


Otherwise 


(2.58) 


V = -2^  . 
E 


(2.59) 


The  above  equation  is  then  numerically  integrated  to  yield  the  recoil  dis- 
tribution. The  results  of  the  above  approach  for  the  case  of  150  ke'/  As^ 
implanted  through  430  X of  Si02  [7]  are  shown  in  Fig.  2.5.  While  the  above 
approach  has  included  a considerable  number  of  approximations,  it  has  been 
shown  to  be  accurate  within  a factor  of  2 for  the  implantation  of  oxygen 
from  thin  Si02  layers  on  Si  bombarded  by  Kr  ions,  and  is  expected  to  give 
at  least  an  order  of  magnitude  estimate  for  the  case  considered  here. 

As  a final  note,  it  is  worthwhile  to  emphasize  the  difficulties 
in  trying  to  develop  theoretical  estimates  of  the  processes  that  occur  in 
implantation  through  a multilayer  target.  We  first  recall  that  range  distri- 
bution calculations  for  single-element  amorphous  targets  produce  no  less 
chan  a 10  - 20%  deviation  from  measured  values,  due  to  such  problems  as 
shell  effects  in  electronic  excitation,  electronic  excitation  during  nuclear 
collisions,  loss  of  ions  through  back  .scattering,  erosion  of  the  target 
surfaces  due  to  sputtering,  and  difficulties  in  developing  accurate  solutions 
to  the  range  equation  (Eq . (2.42))  [9].  If  we  now  consider  the  production  of 
recoils  from  a Si02  layer  of  finite  thickness,  we  an.  io, ed  with  additional 
difficulties.  Theoretical  estimates  of  the  stopping  power  of  SiO.,  give  results 
significantly  in  error  from  experimentally  determined  values  [32].  Measurementi 


of  Si02  films  on  Si  surfaces  indicate  that  nearly  10  oxygen  atoms  are  sputtered 
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Figure  2.5. 

Calculated  depth  distribution  of  oxygen  recoils  for 
a lO^^cm  ^ 150  keV  As  implant  through  430  X of  SiO^ 

(after  Moline  et  £l.  f7]). 
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away  for  each  oxygen  atom  recoil-implanted  for  a 200  keV  As  implant  [5]; 

thus,  a strict  treatment  of  implantation  through  a Si02  layer  would  have  to 

treat  the  target  surface  as  a boundary  chat  moves  as  a function  of  ion  dose. 

This  would  lead  to  a further  spread  in  the  recoil  distribution  beyond  the 

interface  and  would  also  increase  the  recoil  yield  wim  dose,  due  to  the 

fact  that  recoils  from  the  surface  would  have  Co  travel  a shorter  distance 

to  the  substrate  as  the  surface  erodes.  For  example,  analysis  of  an  oxide 

16  "*2 

layer  before  and  after  a 1 x 10  cm  200  keV  As  implantation  yielded  a 
decrease  of  36  A in  an  Si02  layer  originally  176  X thick  [5J  , a 207„  change. 

In  addition  to  affecting  the  recoil  yield,  oxide  erosion  also 
affects  the  shape  of  the  As  distribution  inside  the  Si  layer,  for  exactly 
the  same  reasons.  One  further  consideration  also  affects  the  As  distribution 
inside  the  Si  for  a through-oxide  implant.  Since  has  a greater  stopping 

power  than  Si,  the  distribution  of  As  ions  that  penetrate  the  SiO^-Si  inter- 
face is  peaked  at  a shallower  range  than  for  an  equivalent  implantation  in 
a homogeneous  Si  sample.  This  would  imply  a smaller  penetration  depth 
and  thus  higher  peak  concentration  and  a steeper  fall-off  for  the  As  profile. 
Finally,  the  variation  in  stopping  powers  also  affects  the  mathematical 
treatment  of  implants  through  multiple  layers.  Again,  since  SiO^  has  a 
greater  stopping  power,  ions  that  would  come  to  rest  just  beyond  the  interface, 
based  on  SiO^  range  statistics,  in  fact  travel  further  in  the  Si  due  to  being 
in  a material  with  a lower  stopping  power.  Neglecting  ion  backscattering , 
this  would  yield  a distribution  that  is  discontinuous  at  the  interface  [33]. 
Such  a discontinuity  has  not  been  experimentally  observed. 
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The  theory  of  implantations  into  layered  targets  is  still  in  an 
early  state  of  development,  and  therefore,  experimental  examination  of  a given 
set  of  conditions  provides  the  best  analysis.  However,  as  more  data  becomes 
generated,  we  can  look  forward  to  greater  refinements  in  the  theory. 


j 


3.  niOTOLUMINESCENCK 


Eow-temperature  photo  luminescence  has  been  shown  to  be  an  extremely 
sensitive  monitor  of  radiation  damage  in  silicon  [34-37],  including  damage 
resulting  from  ion  implantation  [38,39].  Luminescence  spectra  from  radiation- 
damaged  silicon  have  been  studied  to  relate  trace  impurities  growth 
techniques  and  the  presence  of  chemical  dopant  impurities  to  the  formation 
and  annealing  of  defects  in  silicon.  For  these  reasons,  it  is  natural  Co 
apply  this  technique  to  the  study  of  the  residual  damage  seen  following 
arsenic  implantations  into  silicon  through  Si02 . The  advantage  of  photo- 
luminescence as  a diagnostic  technique  is  the  sensitivity  of  radiative 
recombination  to  the  presence  of  defects.  It  is  also  possible  to  correlate 
peaks  in  the  recombination  spectra  to  known  defects  which  have  well- 
characterized  impurity  dependences.  Itiis  chapter  contains  an  introductory 
desc.iption  of  Che  technique  of  low-temperature  photoluminescence,  and  the 
impL.rity  dependence  of  the  major  peaks  seen  in  the  recombination  spectra  of 
implanted  silicon.  This  method  is  then  applied  to  the  specific  problem  of 
through-oxide  arsenic  implants. 

3.1.  Photoluminescence  System 

The  luminescence  system  used  in  this  work  is  schematically  illustrated 
in  Fig.  3.1.  Light  from  a 500  W high-pressure  Hg  lamp  is  chopped,  filtered 
to  remove  its  infrared  components,  and  focused  onto  the  surface  of  a sample 
held  in  a liquid  helium  temperature  cryostat.  This  light,  with  photon 
energies  greater  than  the  Si  bandgap,  is  absorbed  at  the  sample  surface, 
creating  excess  carriers.  Some  of  these  carriers  recombine  radiatively,  and 
the  resulting  luminescence  (lower  than  bandgap  energy)  passes  through  the 
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sample  and  exits  at  the  back  surface.  This  luminescence  is  rollected, 
spectrally  analyzed  in  a .larrell-Ash  0.5m  spectrometer,  and  cletectod  by  a PbS 
detector  (d’‘  = 6 x lo'^'^cm(Hz)  /watt)  using  lock-in  techniques.  Spectra 
are  corrected  for  system  response,  detector  output,  monochromator  calibration 
and  sample  observation  temperature.  Detailed  descriptions  of  the  cryostat 
and  optical  system  are  given  elsewhere  [36]. 

3.2.  Silicon  Luminescence  and  Defect  Identification 

Since  silicon  is  an  indirect  semiconductor,  luminescence  at  the 
bandgap  energy  is  not  observed,  and  instead  excitons  recombine  with  the  aid 
of  momentum-conserving  phonons.  This  yields  the  spectrum  shown  in  Fig.  3.2. 
Because  the  TO  phonon  has  the  largest  density  of  states  [40] , this  phonon- 
assisted  band  yields  the  highest  luminescent  intensity.  The  intrinsic 
luminescence  spectrum  shown  in  Fig.  3.2  is  suppressed  by  particle  irradiation, 
and  defect-related  luminescence  is  observed  instead  (Fig.  3.3).  This  lumines- 
cence is  characterized  by  dominant  zero-phonon  peaks  due  to  electronic 
transitions  at  defects  and  by  associated  lower  energy  bands  due  to  phonon 
participation  in  the  transitions. 

The  three  dominant  zero-phonon  peaks  seen  from  ion-implanted  silicon 
are  those  at  1.018  eV,  0.970  eV,  and  0.790  eV.  Hie  1.018  eV  lire  has  been 
identified  [39]  with  a five-vacancy  complex  modelled  by  EPR  [41].  The  0.970  eV 
peak  was  until  recently  identified  with  the  divacancy  [36],  and  the  0.790  eV 
line  was  attributed  to  the  SiG  (15)  (or  K)  center  [36].  However,  the  latter 
two  assignments  have  been  disputed  on  the  basis  of  low-temperature  irradiation 
studies  [37]  and  the  effects  of  A1  and  Ga  impurities  on  these  centers  [42]. 
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Figure  3.2.  Free  exciton  recombination  luminescence  from  unirradiated, 
high  resistivity  Si  showing  bandgap  energy  Eq  and  energy 
of  bandgap  minus  exciton  binding  energy  (after  Haynes  [40]) 
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Figure  3.3.  Recombination  luminescence  from  p-type  pulled  Si  implanted 
with  120  keV  b"^  to  a fluence  of  10^^  ions/cm^  following 
room  temperature  annealing.  The  peak  near  1.10  eV  is  due 
to  free  exciton  recombination,  while  the  peaks  at  1.018  eV  (II), 
0.970  eV  (G),  and  0.790  eV  (C)  are  due  to  irradiation  induced 
defects  (after  Noonan  et  al.  [38]). 


These  recent  studies  suggest  identification  of  the  0.970  eV  line  with  the 

split  illl'  carbon  interstitia Icy  defect,  and  the  0.790  eV  line  with  a 

carbon-modified  A center  (0-vacancy  complex).  Pliotoluminescence  studies 

of  carbon-  or  oxygen-implanted  silicon  were  performed  [43],  and  the  latter 

experiment  suggested  that  impurity  dependence  of  the  0.970  and  0.790  eV 

peaks  could  only  be  determined  by  examining  the  effects  of  incorporated 

impurities  on  the  photoluminescence  spectra  following  irradiation.  For  this 

reason,  (lll'>  n-type,  121-83  0 -cm  silicon  wafers  were  implanted  with 

1 X 10 ^^cm  ^ 150  keV  0 , or  1.1  x lO^^cm  ^ 111  keV  C^,  or  both.  These  doses 

and  energies  were  chosen  to  provide  nearly  identical  impurity  depth  distri- 

18  “3 

butions  in  the  substrates,  with  peak  concentrations  of  5 x 10  cm  [19]. 
These  wafers,  along  with  an  unimplanted  control  wafer,  were  annealed  at  900°C 
under  flowing  argon  gas  to  remove  the  primary  implantation  damage  and  then 
were  irradiated  with  10^^  cm  ^ 2.5  MeV  electrons  at  0°C.  These  wafers  were 
examined  as  irradiated  and  also  after  a 200*^0  anneal,  which  increases  the 
intensity  of  the  peaks.  The  resultant  intensities  of  the  three  major  peaks 
are  summarized  in  Table  I.  From  these  results,  it  can  be  seen  that  the 
0.970  eV  peak  dramatically  increases  in  intensity  for  those  samples  which 
had  been  doped  with  carbon  prior  to  the  irradiation.  After  room  temperature 
annealing,  the  0.970  eV  line  appears  at  greater  intensity  in  the  carbon-  and 
oxygen-doped  substrates  than  in  the  substrate  doped  with  carbon  only;  however, 

following  a 200°C  anneal,  the  intensity  of  this  peak  increases  sharply  for 

« 

the  oxygen-lean  sample.  This  behavior  is  identical  to  that  seen  from  electron 
irradiated  pulled  (oxygen-rich)  and  float  zone  (oxygen-lean)  silicon  [44], 
and  demonstrates  the  interaction  between  carbon  and  oxygen  in  the  recovery  of 
radiation  damaged  silicon  [45].  The  carbon  dependence  of  this  peak  is 
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TABLE  I 


Relative  intensities  of  the  major  luminescence  peaks  from 

ion  implantation-doped,  e lec tron-i rradia ted  float  zone 

(n-type,  121-83  Q -cm)  silicon.  All  implantations  were 

performed  at  7°  off  the  \lll)  axis  at  room  temperature, 

at  energies  of  111  keV  and  150  keV , to  fluences  of  1.1 
14  14  -2 

X 10  and  10  cm  for  carbon  and  oxygen,  respectively. 

The  electron  irradiations  were  performed  at  0°C  at  2.5  MeV 
17  -2 

to  a fluence  of  10  cm  into  material  previously  implanted 
as  indicated  in  the  dopant  column,  followed  by  900°C  anneals. 


Room  Temperature  Anneal  200°C  Anneal 


Dopant 

0.790  eV 

0.970  eV 

1.018  eV 

0.790  eV 

0.970  eV 

1.018 

C,0 

185 

407 

>1. 

332 

383 

* 

C 

13 

330 

* 

92 

580 

* 

0 

110 

75 

* 

187 

119 

17 

29 

49 

■k 

44 

88 

85 

I 


35 


consistent  with  the  previously  suggested  identification  of  this  center  with  the 
split  (ill'  carbon  intersti tia Icy  (Si-G-11)  center  |37J.  Together  with  the 
results  of  the  low  temperature  irradiations  [37]  and  the  effects  of  Ga 
impurities  on  the  spectra  from  irradiated  silicon  [42],  this  work  further 
reinforces  the  proposed  defect  identification. 

The  samples  which  had  been  intentionally  doped  with  oxygen  showed 
very  strong  0.790  eV  signals  compared  with  those  which  had  not  been  oxygen- 
implanted.  This  behavior  is  consistent  with  previous  studies  of  irradiated 
oxygen-doped  silicon  [36].  While  the  oxygen  dependence  of  the  0.790  eV  peak 
has  been  previously  identified,  the  carbon  dependence  of  this  peak  has  not 
been  determined  prior  to  the  work  reported  here.  Table  I indicates  that 
the  0.790  eV  luminescence  is  higher  in  the  samples  doped  with  both  carbon  and 
oxygen  compared  with  those  doped  with  oxygen  alone,  consistent  with  the 
hypothesis  that  the  0.790  eV  luminescence  arises  from  the  carbon-modified  A 
center.  The  weak  presence  of  the  0.790  eV  peak  in  those  samples  which  had 
not  been  oxygen-implanted  appears  to  be  the  result  of  oxygen  coruamination 
from  the  high  temperature  annealing  prior  to  the  electron  irradiation. 

The  growth  of  the  1.018  eV  peak  is  consisuent  with  it.*-  identification 
with  the  five  vacancy  complex.  This  peak  is  seen  strongly  for  ion-implanted 
silicon  after  room  temperature  irradiations,  but  has  never  been  seen  directly 
after  electron  irradiations  [43].  That  this  center  is  seen  in  these  electron- 
irradiated  samples  is  an  effect  of  the  interaction  of  impurities  with 
irradiation-induced  defects  during  the  annealing  process.  Both  carbon  and 
oxygen  are  efficient  vacancy  traps  [9.J.  Simpler  defeeth,  such  as  the  E (.donor- 
vacancy)  center,  break  up  by  200°C,  liberating  vacancies  which  can  in  turn 
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hecome  trapped  by  carbon  or  oxygen  impurities.  One  would  then  expect  to  see 
the  highest  intensity  from  a vacancy- re  la  ted  center  in  those  samples  which 
have  the  lowest  carbon  and  oxygen  concentrations.  Thus,  we  see  in  Table  I 
that  the  1.018  e'.’  peak  is  greatest  in  intensity  at  200”c  annealing  for  the 
control  sample,  and  next  largest  in  intensity  for  the  sample  doped  with 
oxygen.  This  behavior  correlates  well  with  the  already  strong  case  for  the 
identification  of  the  1.018  eV  peak  with  the  five  vacancy  complex  [39]. 

In  addition  to  the  peaks  mentioned  above,  new  peaks  were  seen  in 
the  photoluminescence  spectra  from  the  irradiated  samples.  Unfortunately, 
the  intensities  of  these  peaks  were  too  small  to  allow  detailed  investigation. 
These  samples,  with  their  dopant  and  annealing  dependencies,  are  summarized 
in  Tab le  II . 

3.3.  Oxygen  Recoil  Detection 

With  further  evidence  establishing  the  0.790  eV  peak  as  arising 
from  an  oxygen-dependent  center  and  the  1.018  eV  peak  as  arising  from  the  five 
vacancy  complex,  we  are  now  in  a position  to  apply  this  knowledge  to  the 
study  of  through-oxide  arsenic  implantations  into  silicon.  In  this  study, 
bare  (no  grown  oxide)  float  zone  (oxygen- lean)  (111''  n-  and  p-type  wafers 
1 123-81  n cm)  were  Implanted  with  As  ions  to  isolate  the  damage  due  to  As 
implantation  alone.  To  determine  the  effects  of  possible  oxygen  recoils, 
similar  float  zone  crystals  were  implanted  through  a thin  (~  500  layer 
of  SiO^  grown  by  dry  thermal  oxidation.  As  a control  to  test  the  effects  of 
oxygen  in-diffusion,  an  identical  float  zone  sample  was  oxidized  under  the 
same  conditions  as  the  above  wafers,  and  then  the  Si02  was  removed  by  chemical 
etching  prior  to  implantation. 
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TABLE  II 

Additional  zero-phonon  peaks  seen  from  electron-irradiated 
silicon  samples,  listing  the  dopants  added,  if  any,  and  the 
anneals  (room  temperature,  200°C,  or  both)  for  which  they 
were  observed. 

Peak  Energy  (eV)  Samples  and  Anneals 

0.950  C and  0 (RT,  200°C) 

0.912  C and  0 (RT,  200°C) 

0.900  Undoped  (200°C) , 0 (RT , 200°C) 

0.898  C and  0 (RT , 200°C) , (RT,  200°C) 

0.890  C and  0 (200°C) , 0 (RT,  200°C),  undoped  (200°C) 

0.811  0 (RT) 
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Arsenic  implantations  of  10  cm  were  then  performed  at  room 


temperature  at  200  keV  with  typical  beam  current  densities  of  0.05  y A/cm 


Wafers  were  oriented  7 from  the  beam  direction  to  minimize  channeling 


effects.  Anneals  below  300  C were  performed  in  a controlled  temperature  oil 


bath,  whereas  anneals  above  300  C were  performed  in  a tubular  funvico  under 


flowing  argon  gas.  All  anneals  in  ttiis  study  were  of  30  minutes  duration 


We  examine  first  the  luminescence  p-.  aks  common  to  all  spectra 


seen  from  the  arsenic- implanted  samples.  These  peaks  occur  at  1.018  eV  (11), 


1.034  (14),  and  1.039  eV  (13).  Of  all  the  peaks  examined  in  this  study,  the 


peak  at  1.018  eV  produced  the  greatest  intensity.  This  peak  grew  in 


following  the  150°C  anneal,  reached  an  intensity  maximum  following  250°C 


annealing,  and  then  disappeared  following  350  C anneals.  The  1.018  eV  peak 


appears  at  such  great  intensity  that  for  the  200°C  and  250°C  annealing 


it  completely  dominates  the  recombination  spectrum,  appearing  at  approxi- 


mately 40  times  the  intensity  seen  from  the  electron- irradiated  samples 


(Table  I).  The  intensity  of  this  peak  is  such  that  we  are  able  to  identify 


its  replication  by  the  18  meV  TA,  the  44  meV  LA,  the  53  meV  LO,  and  the 


59  meV  TO  phonons  (Fig.  3.4).  The  14  peak  1.034  eV  appears  as  a higher 


energy  shoulder  on  the  1.018  eV  peak  for  the  200  C anneal,  and  then  vanishes 


for  higher  temperature  annealing.  Peak  13  (1.039  eV)  appears  weakly  for 


the  300°C  anneal,  reaches  a maximum  intensity  for  350°C,  and  then  anneals 


out  following  450  C annealing.  The  annealing  behavior  of  all  three  of 


these  peaks  is  identical  to  that  observed  from  samples  implanted  with  other 


species  [39],  and  indicates  that  these  peaks  are  probably  due  to  defects 


without  impurity  interaction  or  involving  trace  impurities  common  to  silicon 


substrates . 


I 


Relative  Intensity 


Photon  Energy  (eV) 


Figure  3.4.  Photoluminescence  from  float  zone  Si  following  As 

implantation  through  an  Si02  layer.  The  intense  peak 
at  1.018  eV  is  due  to  a five  vacancy  cluster.  The 
weak  peak  near  0.79  eV  is  due  to  an  oxygen  dependent 
defect.  The  0.790  eV  peak  is  absent  in  the  sample 
implanted  after  removing  the  SiO,  layer  (broken  line). 
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The  most  Interesting  results  of  tliis  study  come  from  examination 
of  the  0.79  eV  peak.  This  peak  is  seen  from  all  the  through-oxide  implanted 
float  zone  (oxygen- lean)  samples  as  well  as  from  the  bare  surface  implanted, 
pulled  (oxygen-rich)  substrates,  but  not  from  the  bare  surface  implanted, 
float  zone  samples  Furthermore,  the  control  float  zone  sample,  from 
which  the  SiO^  was  stripped  before  implantation,  exhibited  no  oxygen- re  la  ted 
luminescence  (Fig.  3-4).  The  presence  of  this  oxygen-dependent  damage 
ci-nter  for  the  float  zone  samples  which  had  been  implanted  through  a Si02 
layer,  and  not  for  the  bare  surface  implantations  nor  for  the  control 
sample,  siiows  conclusively  the  arsenic  implantations  through  silicon 
dioxide  result  in  oxygen  recoils  which  enter  the  silicon  substrate.  These 
oxygen  atoms  participate  in  irradiat ion- induced  defect  complexes,  including 
that  related  to  the  0.79  eV  luminescence.  It  is  seen  that  oxygen  in-diffusion 
during  thermal  growth  of  the  SiO^  is  insufficient  to  produce  detectable 
luminescence  following  arsenic  implantations. 

The  oxygen-dependent  damage  center  causing  the  0.790  eV  lumines- 
cence anneals  out  at  350°C,  so  that  this  center  is  not  itself  responsible 
for  the  residual  damage  seen  following  through-oxide  implantations.  However, 
the  presence  of  the  0.790  eV  luminescence  provides  a direct  detection  of 
the  oxygen  recoil,  and  allows  us  to  estimate  an  average  oxygen  concentration. 


The  intensity  of  the  0.790  eV  peak  for  the  through-oxide  implanted  samples 


was  half  that  of  the  arsenic- implanted  pulled  silicon  sample.  Since  the 

18  - 3 

oxygen  content  of  pulled  silicon  is  typically  1 x 10  cm  [^7  “ , by  assuming 


i 


the  production  rate  of  this  defect  is  approximately  the  same  for  recoil- 
implanted  oxygen  as  for  grown  silicon  --  an  assumption  which  may  not  be 
strictly  justified  [43'  --  we  estimate  the  "average"  oxygen  concentration 


41 


4 


17  - 3 

ill  till'  ri'coi  1 - imp  lant  I'd  layer  to  be  on  tlie  order  ol  5 x 10  cm  . This 
I'Stimate  agrees  well  with  the  calculations  of  Moline  jet  (Figure  5)  if 

their  concentration  is  scaled  down  to  represent  an  implantation  at  thi- 
appropriate  fluence. 

The  samples  were  annealed  beyond  400°C  to  compare  the  spectra 
from  the  through-oxide  implanted  samples  to  those  from  the  bare  surface 
implants.  The  luminescence  disappeared  by  500°C  and  did  not  return  even 
after  1000°C  annealing,  except  for  weak  luminescence  from  the  unimplanted 
substrate.  For  this  reason,  it  is  necessary  to  resort  to  other  techniques 
to  study  the  residual  damage  resulting  from  the  through-oxide  implants. 

The  next  chapters  will  examine  and  compare  diffusion,  electrical  activation, 
and  carrier  mobility  profiles  for  through-oxide  and  bare  surface  As 
implantations . 

It  should  be  noted  in  closing  that  MOS  techniques  have  also  been 
used  to  detect  the  oxygen  recoil  [47].  This  method  relies  on  the  presence 
of  surface  states  with  well  characterized  energy  levels  in  the  silicon  bandgap 
to  provide  defect  identification.  It  is  reported  that  the  density  of  these 
surface  states  can  be  calibrated  to  yield  a direct  measure  of  oxygen  recoil 
densities  (47),  the  results  of  which  are  again  in  order  of  magnitude  agreement 
with  the  calculations  of  Moline  [7]. 
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4 . ELKMKNTAL  PROFILING 

Not  all  implanted  ions  become  electrically  active  after  annealing; 
the  distribution  of  electrically  active  impurities  differs  from  the  total 
impurity  distribution  due  to  the  presence  of  defects.  Thus,  the  only 
accurate  method  of  evaluating  the  activation  of  implanted  ions  is  to 
compare  their  elemental  profiles  with  the  profile  of  electrically  active 
impurities.  This  chapter  will  be  concerned  with  the  measurement  of  elemental 
depth  profiles  for  implanted  arsenic  and  recoil  implanted  oxygen  and  carbon. 

The  next  chapter  will  deal  with  measurement  of  the  electrical  properties 
of  the  implanted  layers  and  comparison  of  the  carrier  and  mobility  profiles 
with  the  elemental  impurity  profiles. 

4.1.  Sample  Preparation 

The  wafers  used  for  the  determination  of  electrical  and  elemental 
profiling  were  all  cut  from  (lOO^  - oriented  boron  doped  single  crystal 
silicon  boules  and  all  had  commercially  polished  surfaces.  One  set  of  5,2 
to  7 O'cm  Ventron  substrates  was  implanted  with  200  keV  As"*"  to  a fluence 
of  1 X iO^^cm  ^ cither  into  the  bare  surface  or  through  approximately  840  R 
of  thermally  grown  Si02  (1100°C,  30  min  dry  oxygen).  The  other  set  of  wafers 

used  in  this  study  was  Texas  Instruments  9-14  fl-cm  pulled  silicon.  These  wafers 

+ ! 5 -2 

wore  Iu.j  lai>Led  with  150  ko'.‘  As  to  a fluence  of  5 x 10  cm  , <.  tuner  into 
the  bare  surface  or  through  approximately  150  R of  silicon  dioxide 
(900*^0,  40  min  dry  oxygen).  These  samples  were  examined  as  implanted  and 
also  after  selected  annealing  stages.  All  anneals  in  this  study  were  performed 
in  a quartz-lined  tubular  furnace  under  flowing  argon  gas.  line  samples 
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implanted  to  5 x 10  cm  “ As  were  also  examined  following  redistribution 
by  thermal  oxidation.  Both  the  bare  surface  and  the  through  oxide  implanted 
samples  were  subjected  to  either  5 minutes  anneal  ir.  dry  nitrogen,  2 hours 
wet  oxidation,  and  5 minutes  anneal  in  dry  nitrogen  performed  sequentially 
at  900  C,  or  5 minutes  annea  1 in  dry  nitrogen,  1 hour  wet  oxidation  and 
5 minutes  anneal  in  dry  nitrogen  performed  sequentially  at  1100°C.  Oxidc- 
thickness  values  quoted  in  this  study  were  measured  by  e llipsomet ry . Table  111 
summarizes  the  samples  used  for  profile  determination. 

4.2.  Techniques  for  Elemental  Profiling 

Three  techniques  were  available  in  this  study  for  the  determination 
of  elemental  profiles;  Auger  electron  spectroscopy  (AES)  [48],  glow  discharge 
optical  spectroscopy  (GDOS)  [49],  and  secondary  ion  mass  spectroscopy  (SIMS) 
[50].  These  techniques  all  employ  sputtering  for  sample  sectioning,  but 
differ  in  the  ways  in  which  the  profile  information  is  obtained,  llie  mechanisms 
basic  to  these  techniques  will  be  briefly  summarized  below. 

4.2.1.  Auger  Electron  Spectroscopy 

In  this  technique,  the  sample  to  be  examined  is  bombarded  in  high 
vacuum  (<"  10  torr)  by  a primary  electron  beam  in  the  energy  range  of 

1-10  keV . This  primary  beam  causes  a variety  of  secondary  processes  to  occur, 
the  most  important  of  which,  for  this  technique,  is  the  creation  of  an 
atomic  inner  shell  vacancy.  ITie  inner  shell  vacancy  is  filled  by  an  outer 
shell  electron,  and  this  excess  energy  is  carried  away  either  by  a photon 
or  by  an  electron.  If  an  electron  is  ejected,  the  radiationless  transfer  of 
energy  is  known  as  the  Auger  effect  and  the  emitted  electron  as  the.  Auger 
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Summary  of  samples  used  in  this  study  to  determine  elemental 
and  electrically  active  profiles.  All  substrates  used  were 
( 100''  oriented  p-type  (boron-doped) . Values  quoted  for  oxide 
thickness  were  determined  by  e 1 li psome try . 


Pre-Implant 

Sample 

Sample 

Surface 

As  Energy 

Post-Implant 

Tvpe 

Resistivity 

Condition 

and  Dose 

Oxidation 

1 

9 - 14  n-cm 

~ 150 S SiO^ 

C -a15  -2 

5xi0  cm 

150 

keV 

None 

2 

9-14  Q-cm 

~ 150 X SiO. 

. ,pl5  -2 
5x10  cm 

150 

keV 

2 hr  <2 

<3  900°C 

1 

3 

9-14  n-cm 

~ 150  SiO, 

5xl0‘^cm  ^ 

150 

keV 

1 hr  wet  0, 
(?  1100°C 

4 

9-14  fi-cra 

Bare 

^ -nl5  -2 
5xi0  cm 

150 

keV 

2 hr  wc  - 02 
'a  900°C 

5 

9-14  Q-cm 

Bare 

,^15  -2 

5x10  cm 

150 

keV 

None 

6 

9 - lit  n-cm 

Bare 

- -mis  -2 
3x10  cm 

150 

keV 

1 hr  w C'2 

(a900°C 

13 

5.2-7  n-cm 

~ 840^  SiO^ 

ixlO^^cm  ^ 

200 

keV 

None 

14 

5.2-7  n-cm 

Bare 

...15  -2 

ixlO  cm 

200 

keV 

None 

45 


electron  1 52 1 . I'ht  energy  of  this  ejected  electron  is  determined  by  the 
energy  levels  of  the  atom  before  and  after  ejection  and  hence  is  characteristic 
of  the  emitting  atom.  If  the  emitting  atoms  are  close  enough  to  the  sample 
surface  (3-20,^  f51]),  the  emitted  electrons  will  suffer  no  inelastic  collisions 
and  tlius  retain  their  characteristic  energy  Electrons  emitted  from  greater 
depths  will  suffer  inelastic  collisions  and  thus  lost:  their  characteristic 
energies.  These  ine lastica 1 ly  scattered  Auger  electrons  combine  with 
secondary  electrons  created  by  other  mechanisms  to  form  a broad  continuum  of 
electrons  upon  which  the  Auger  peaks  are  superimposed.  In  this  technique, 
the  voltage  of  the  electron  analyzer  is  modulated  by  a few  volts  and  lock-in 
detection  is  used  to  electronically  differentiate  the  low  level  Auger  peaks 
from  the  large  continuous  background  [53].  Since  only  the  Auger  electrons 
emitted  very  close  to  the  sample  surface  retain  their  characteristic  energies, 
Auger  spectroscopy  is  essentially  a surface  analytical  techiiique.  Depth 
profiles  are  obtained  by  gradually  sputtering  the  sample  surface  to  expose 
successively  deeper  regions  to  Auger  analysis. 

4.2.2.  Glow  Discharge  Optical  Spectroscopy 

In  the  GDOS  technique,  the  sample  to  be  examined  is  placed  on  the 
cathode  of  a dc  glow  discharge  operated  in  the  "abnormal  glow"  mode  [54]. 

Atoms  from  the  sample  to  be  examined  are  ejected  irtc  the  discharge  by  the 
sputtering  gas  (typically  argon  or  krypton).  Most  of  tnese  sputtered  atoms 
leave  the  sample  surface  as  neutrals  [55].  Some  of  these  atoms  are  collisiona lly 
excited  by  secondary  electrons  escaping  the  cathode,  and  emit  characteristic 
atomic  luminescence.  Ttiis  luminescence  is  spectrally  analyzed  for  the  emission 
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characteristic  of  the  element  of  interest,  and  the  intensity  of  this  lumines- 
cence is  monitored  with  time.  llie  GDOS  process  assumes  that  the  emission 
intensity  is  proportional  to  the  concentration  of  the  element  of  interest 
in  the  substrate,  llius , emission  intensity  can  be  calibrated  by  use  of  suitable 
standards  to  yield  impurity  concentration,  and  the  sputtering  rate  can  be 
independently  determined  to  yield  a depth  calibration.  It  is  important  to 
note  that  the  element  to  be  profiled  is  examined  in  a region  outside  the 
sample  surface,  thus  reducing  the  influence  of  the  substrate  on  the  profile 
measurement  to  that  due  to  the  mechanisms  of  the  atomic  ejection  during 
sputtering . 

4.2.3.  Secondary  Ion  Mass  Spectroscopy 

The  SIMS  technique  employs  a focused  primary  ion  beam  to  bombard 
the  sample  to  be  examined.  Unlike  GDOS,  however,  SIMS  analyzes  the  small 
fraction  of  sample  atoms  that  leave  the  surface  as  ions.  These  ions  are 
extracted  from  the  sputtered  area  and  mass  analyzed  to  provide  the  profile 
information.  As  in  GDOS,  however,  the  atomic  concentration  is  produced  as  a 
direct  result  of  the  sputtering  process;  thus  both  GDOS  and  SIMS  are  inherently 
depth  profiling  techniques. 

As  with  Auger  spectroscopy,  the  escape  depth  for  the  analyzed 
particles,  here  the  secondary  ions,  is  10-20  R [51];  however,  the  interpreta- 
tion of  SIMS  data  is  significantly  more  complicated  than  is  AES  data.  The 
reason  for  this  complication  lies  in  the  various  natures  of  the  two  techniques. 
Auger  electrons  are  produced  by  atomic  inner  shell  vacancies;  the  energies  of 
these  core  levels  are  barely  affected  by  chemical  binding.  Ionization,  however. 
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occurs  at  the  outer  shell  of  electrons  --  where  all  the  chemical  binding 
occurs.  Thus  Che  secondary  ion  yield  can  vary  over  several  orders  of 
magnitude  for  different  elements  in  the  same  matrix  or  for  the  same  element 
in  different  matrices  [50].  Typically,  the  secondary  ion  yield  of  an 
element  is  much  higher  when  it  is  sputtered  from  a bound  state  than  is  the 
case  for  a metallic  state  [50].  Thus,  the  presence  of  a reactive  species 
in  Che  sample  enhances  the  secondary  ion  yield  [50].  SIMS  profiling 
typically  employs  either  oxygen  or  cesium  as  the  primary  ion  beam  to  stabilize 
the  secondary  ion  yield  at  high  levels.  From  basic  chemistry,  we  recall 
that  elements  in  the  left  of  the  periodic  table  (such  as  boron)  tend  to  lose 
electrons  to  form  complete  outer  shells;  these  elements  thus  prefer  to 
form  positive  ions,  and  therefore  oxygen  (electron-accepting)  primary  beams 
are  used  to  enhance  secondary  ion  yields  for  these  elements.  On  the  other 
hand,  elements  from  the  right  of  the  periodic  table,  such  as  arsenic,  attract 
electrons  to  complete  their  outer  shells  and  form  negative  ions;  thus  cesium 
(electron-donating)  primary  beams  are  used  to  enhance  the  secondary  yields 
of  these  elements  [56].  Nevertheless,  by  examining  know.T  concentrations  of  a 
particular  element  in  a given  matrix,  calibration  of  sample  concentration  can 
be  obtained,  and  as  in  GDOS  the  sputtering  rate  can  be  independently  measured 
to  provide  a direct  depth  calibration. 

4.3.  Results 

4.3.1.  Auger  Electron  Spectroscopy 

The  Auger  profiles  of  carbon,  oxygen,  and  arsenic  for  the 
! 5 -2 

5 X 10  cm  150  ke . As-implanted  samples  are  shown  in  Fig.  4.1.  Ihis  data 
v;as  obtained  on  a PHI  system  using  a 5 keV  coaxial  electron  gun. 
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Figure  4.1.  AES  profiles  for  oxygen,  carbon,  arsenic  and  silicon, 

15  -2 

for  unannealed  5 x 10  cm  , 150  keV  As  implants  into 
bare  surface  (top)  and  through  Si02  (bottom).  Notice 
lint  r intensity  scales. 
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It  can  be  seen  that  even  for  the  peak  As  concentration  of  7 x 10  cm 
expected  for  the  bare  surface  implant  from  LSS  theory  [19],  the  signal  appears 
very  noisy.  Auger  spectroscopy  has  very  pocr  sensitivity  (minimum  detectable 
concentration  0.17o  of  matrix  for  elements  with  high  atomic  numbers  [51]), 
although  better  sensitivity  can  be  obtained  for  elements  of  lower  atomic 
number,  such  as  carbon  and  oxygen.  Unfortunately,  however,  surface  contamina- 
tion by  these  common  trace  elements  is  unavoidable,  so  that  Auger  spectra 
always  show  strong  initial  carbon  and  oxygen  signals  even  when  the  native 
oxide  is  chemically  removed  immediately  prior  to  insertion  in  the  Auger  system. 
As  the  sample  surface  is  sputtered  to  develop  the  depth  profiles,  a fraction 
of  this  surface  layer  is  recoil-implanted  into  the  substrate,  thus  producing 
a sloping  profile  from  the  surface  that  obscures  the  near  surface  carbon  and 
' oxygen  profiles. 

The  only  samples  which  showed  sufficient  arsenic  concentrations  to 
be  easily  measured  were  the  samples  wet  oxidized  at  900°C  after  the  implant. 

For  these  samples,  as  the  oxide  is  grown  through  the  implanted  layers,  the 
arsenic  accumulates  at  the  interface.  Apparently,  at  900^C  the  diffusion  of 
arsenic  away  from  the  interface  is  slower  than  the  advance  of  the  interface 
with  the  wet  oxide  growth,  and  thus  the  entire  implanted  arsenic  dose  Is 

contained  withir,  a region  only  a few  hundred  Angstroms  from  the  interface. 

15  -2  + 

Fig.  4.2  shows  the  Auger  results  for  the  5 x 10  ' cm  150  keV  As  bare  surface 
I.  implant,  after  subsequent  wet  (..xidatic:.  at  900*^0.  Of  special  interest  is 

the  sample  from  which  the  Si02  had  been  chemically  removed  before  examination. 
Of  the  three  elemental  profiling  techniques  used  here.  Auger  spectroscopy  is 
best  suited  for  samples  where  t'--  i.r.purity  distribution  is  peaked  at  th. 
surf.ice,  since  it  does  not  require  the  establishment  of  an  equiUbrinm  in  the 
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Figure  4.2.  AES  results  for  5 x 10  cm  , 150  keV  As  implant 

following  2 hour  wet  oxidation  at  900°C.  For  the  top 
profile,  the  SiO^  was  chemically  removed  prior  to 
examination;  for  the  bottom,  the  profiles  were  taken  by 
sputtering  through  the  SiO^  layer.  Apparent  in  the 
lower  profiles  are  the  loss  of  intensity  due  to  oxide 
charging  effects  and  the  decrease  in  depth  resolution 
with  sputtering  time. 
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sputtering  process  before  a reliable  signal  can  be  developed.  Unfortjnately , 
however,  the  Auger  intensity  obtained  from  these  samples  could  not  be  converted 
into  concentra t ion  due  to  lack  of  a standard  sample  with  sufficiently  high 
arsenic  concentration.  Additionally,  Fig.  4.2  graphically  illustrates  a 
disadvantage  common  to  all  sputtering  techniques;  namely,  the  gradual  loss 

of  depth  resolution  with  .sputtering  time  due  to  nonu;;.'.i\ t;n  s pu ..t . ..g  

[57],  In  addition,  the.  sample  examined  through  the  thick  oxide  also  exhibits 
a lower  Auger  yield  due  to  trapped  charges  in  the  oxide  repelling  some  of 
the  emitted  electrons.  Due  to  all  the  problems  listed  above,  intensity  vs. 
time  results  were  not  calibrated  to  yield  concentration  vs.  depth  plots  for 
the  Auger  analysis. 

4. 3 2.  Glow  Discharge  Optical  Spectroscopy 

Due  to  the  presence  of  hydrocarbon  pump  oils  in  the  discharge 

chamber,  the  GDOS  technique  is  unable  to  provide  elemental  profiles  for  carbon 

or  oxygen.  Itie  elemental  As  implants  are  shown  in  Fig.  4.3  for  the  bare 

surface  and  through  oxide  implanted  samples.  The  bare  surface  impla;-t  shows 

a projected  range  of  820  + 2C  S with  a standard  deviation  for  the  Gassian 

approximation  of  335  + 20  X.  Hiese  values  compare  well  with  compiled  values 

for  As  range  stati.stics  of  842  X projected  range  and  292  X standard  deviation 

[19].  It  is  also  seen  from  this  figure  that  the  detection  limit  for  arsenic 

19  -3 

using  this  technique  is  approximately  3 x 10  cm  . In  the  GDOS  technique, 
the  sensitivity  is  limited  by  the  sputtering  rate  [58];  however,  faster 
sputtering  rates  increase  depth  uncertainty  for  shallow  profiles  such  as  these. 
It  was  determined  during  this  study  that  a sputtering  rate  of  50  8/minute 
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Figure  4.3.  Arsenic  profiles  by  GDOS  for  5 x 10  cm  150  keV  As 
implants.  The  solid  line  represents  the  profile 
resulting  from  a bare  surface  implant  while  the  broken 
line  represents  the  profile  resulting  from  implants 
through  an  approximately  150  K thick  Si02  layer. 
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produced  optimal  results. 

When  the  glow  discharge  is  initiated,  there  is  a temporary  period 
of  instability  until  an  equilibrium  sputtering  rate  is  attained.  This 
initial  instability  could  result  in  an  uncertainty  in  depth  calibration 
for  the  GDOS  process  if  the  glow  discharge  were  initiated  on  the  sample 
surface.  To  avoid  these  problems,  a thin  (approximately  200  X)  nickel  layer 
is  evaporated  over  the  sample  surface,  and  the  discharge  is  initiated  on  this 
nickel  layer.  Nickel  was  chosen  for  this  layer  due  to  its  excellent  adhesion 
to  silicon,  its  faster  sputtering  rate  than  the  silicon  substrate  (thus  avoiding 
cone  formation  [57]),  and  its  excellent  mechanical  properties  (i.e.  scratch 
resistance).  The  present  GDOS  system  employs  two  monochromator  and  detection 
systems  operating  simultaneously.  One  system  monitors  the  intensity  of  the 
element  to  be  profiled,  and  the  second  monitors  substrate  (i.e.  silicon) 
intensity.  This  second  system  indicates  the  time  at  which  the  silicon  surface 
is  reached  and  provides  a monitor  of  the  uniformity  of  the  substrate  sputtering 
rate . 

Problems  were  encountered  In  measuring  the  profiles  of  the  other 
arsenic-implanted  samples.  Difficulties  in  depth  calibration  were  encountered, 

+ o 

however, for  the  150  keV  As  implants  annealed  to  1000  C,  and  uniform 
sputtering  rates  could  not  be  obtained  for  the  sample  redistributed  by  thermal 
oxidation.  Unfortunately,  not  enough  material  from  200  keV  As"*"  implants  was 
available  to  permit  GDOS  measurement. 

4.3.3.  Secondary  Ion  Mass  Spectroscopy 

Since  only  a small  amount  of  sample  area  is  consumed  in  the  SIHS  Loch- 
2 

nique  (approximately  0.01  mm  prbfile)  all  the  samples  mentioned  in  Table  ill 


54 


were  profiled  by  the  SIMS  technique.  On  the  SIMS  system  used  in  this  study 

(an  AEI  model  IM  20  modified  by  the  addition  of  a cesium  ion  source)  it  was 

impossible  to  simultaneously  profile  more  than  one  element;  thus  profiles  of 
different  elements  in  the  same  sample  are  obtained  by  consecutive  runs  from 
areas  immediately  adjacent  on  the  sample  examined.  At  least  two  craters 
are  necessary  for  every  sample  examined--one  to  measure  the  arsenic  intensity 

and  the  other  to  measure  silicon  intensity,  lliis  procedure  was  required, 

since  ion  yield  varies  from  one  position  in  the  sample  holder  to  the  other, 
due  to  geometric  factors.  These  factors  were  fairly  constant  for  each 
particular  sample  position  in  the  holder,  and  thus  could  be  cancelled  by 
normalizing  the  observed  arsenic  intensity  to  the  intensity  of  the  silicon 
substrate . 

The  SIMS  technique  yields  plots  cf  ion  yield  vs.  sputtering  time, 

as  shown  in  Fig.  4.4.  These  plots  were  then  calibrated  to  yield  graphs  of 

concentration  vs.  depth  using  a procedure  described  below.  Arsenic  yield 

normalized  to  the  silicon  yield  was  calibrated  by  assuming  that  the  peak 

15  -2 

value  obtained  for  the  5 x iO  cm  150  keV  bare  surtace  implant  was 
20  -2 

7.0  X 10  cm  (the  value  given  by  LSS  calculations  [19j).  This  value  agrees 

very  closely  with  the  GDOS  result  shown  in  Fig.  4.3.  Depth  scales  were 

calibrated  through  the  use  of  a Sloan  Dektak  to  determine  crater  depths  for 

known  sputtering  times,  and  thus  determine  average  sputtering  rates  for  four 

15  -2  + 

of  the  samples  examined  (5  x 10  cm  150  keV  As  implants  into  bare  surface  or 
through  silicon  dioxide  for  the  600°C  and  1000*^0  anneals).  Since  the  measure- 
ment of  a crater  100  urn  on  a side  with  the  Sloan  Dektak  is  difficult  at  best, 
the  sputtering  rate  was  also  calibrated  by  assuming  that  the  peak  of  the 
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Figure  4.4.  Plot  of  sputtering  As  secondary  ion  yield  (log  scale) 

vs.  sputtering  time  following  600°C  annealing  of 
15  -2  + 

5 X 10  cm  150  keV  As  implant  through  approximately 
150  ^ of  Si02 • This  raw  data  was  converted  to  yield 
arsenic  concentration  profiles  through  the  use  of 
calibration  procedures  discussed  in  the  text. 
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arsenic  distribution  lor  the  150  keV  bare  surface  arsenic  implant  occurred 
at  the  hSS  value  of  842  X [19].  The  difference  bc-tveen  thtse  two  measurement 
techniques  was  less  than  10/,  for  ti’ose  samples  in  wliich  crater  depth  could 
be  unambiguously  determined. 

The  plots  of  arsenic  concentration  vs.  depth  obtained  by  this 
technique  are  shown  in  Figs.  4.5-4.12.  An  arsenic  profile  fr.r  the  600°C 
annealed  sample,  for  which  there  is  negligible  diffusion  from  the  initial 
implanted  profile,  is  shown  in  Fig.  4.5.  Comparing  these  results  with 
the  profiles  obtained  following  lOOO^C  anneal  (Fi  , . ) , we  observe  of  the 

effects  of  concentration  dependent  diffusion  of  arsenic  in  silicon 
[1,59|  for  the  doses  used  here.  Thi=  effect  is  due  to  field-aided  diffusion  [3J. 
Because  of  the  concentration  dependent  diff-sic;, , the  peak  of  the  implanted 
distribution,  having  a higher  arsenic  conceutra ti on , diffuses  more  rapidly  and 
eventually  overtakes  the  tail  of  the  distribution  which  has  a lower  arsenic 
concentration  and  thus  lower  diffusion  coef  f ioie.nt . iTiis  produces  an  abrupt 
arsenic  prefile,  wh.ich  is  extremely  useful  for  emitter  formation  in  bipolar 
transistors  [1].  Tb.is  effect  is  also  seen  ftr  the  samples  impla-ited  through  an 
approximately  150  X thick  SiC2  layer,  where  the  peak  of  the  arsenic  implant  is 
inside  the  silicon  substrate.  T!;e  effe.'.t  is  nut  easily  observable  for  the 
.samples  implanted  through  840  A of  SiO,  , where  the  arsenic  distribution  was 

peaked  near  the  silicon-silicon  dioxide  interfac,t.  In  fact,  for  the  latter  case 

o o 

the  1000  C profile  (Fig.  4.10)  is  s'l’cn  shallower  than  that  of  the  602  i.’ 

annealed  sample  (Fig.  4.9).  This  re.sult  for  implants  through  the  tlilcler 

silicon  dioxide  layer  sugge.sts  that  the  thin  oxide  layer.s  grown  were 


inadvertently  of  non-uniform  thickness.  Subsequent  investigation  by  GDOS , SIMS, 
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Figure  4.5.  Arsenic  profile  by  secondary  ion  mass  spectroscopy 
following  600°C  30  minute  annealing  of  5 x lO^^cm 
150  keV  As  bare  surface  implant. 
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Figure  4.6.  Arsenic  profile  by  secondary  ion  mass  spectroscopy 
following  1000°C  30  minute  annealing  of  5 x lO^^cm 
150  keV  As"^  bare  surface  implant. 
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Figure  4.7.  Arsenic  profile  in  silicon  by  secondary  ion  mass 

spectroscopy  following  600°C  30  minute  annealing  of 
15  -2  + 

5 X 10  cm  150  keV  As  implant  through  an  approximately 

150  R thick  Si02  layer. 
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Figure  4.8.  Arsenic  profile  in  silicon  by  secondary  ion  mass 
spectroscopy  following  1000  C annealing  of  5 x 10 
cm”^  150  keV  As^  implant  through  an  approximately 
150  X thick  SiO.,  layer. 


Figure  4.10.  Arsenic  profile  in  silicon  by  secondary  ion  mass 

spectroscopy  following  1000  C 30  minute  annealing  of 
1 X lO^^cm'^  200  keV  As"*"  implant  through  an 
approximately  840  X thick  SiO^  layer. 
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Figure  4.11.  Arsenic  profile  in  silicon  by  secondary  ion  mass 

15  -2  + 

spectroscopy  for  5 x 10  cm  150  keV  As  implant 
following  redistribution  by  wet  oxidation  for  60 
minutes  at  1100°C.  Profiles  are  identical  for 
implants  into  bare  surface  and  implants  through 
approximately  150  S Si02. 
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Figure  4.12.  Arsenic  profile  in  silicon  by  secondary  ion  mass 

15  -2  + 

spectroscopy  for  5 x 10  cm  150  keV  As  implant 
following  redistribution  by  wet  oxidation  for 
120  minutes  at  900°C.  Profiles  are  identical  for 
implants  into  bare  surface  and  implants  through 
approximately  150  X Si02" 
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and  electrical  measurements  described  in  the  next  chapter  indicated  variations 
in  the  projected  ranges  for  arsenic  implanted  through  the  silicon  dioxide 
layers.  These  variations  were  greater  than  the  uncertainties  of  any  of  the 
depth  profiling  techniques  used.  Such  variation  in  the  arsenic  distributions 
was  not  seen  for  the  bare  surface  implants,  and  thus  could  only  have  arisen 
from  variations  in  the  thickness  of  the  grown  silicon  dioxide  layers,  since 
the  bare  surface  and  through-oxide  implanted  layers  were  implanted 
simultaneously.  Tlie  variation  in  projected  ranges  for  the  th rough -oxide 
implanted  samples  over  various  regions  of  the  implanted  wafer  makes  it 
impossible  to  quantitatively  compare  the  diffusion  of  through-oxide  implanted 
silicon  with  bare  surface  implanted  silicon. 

It  is  interesting  to  observe  the  effects  of  thermal  oxidation  on 
the  arsenic  implanted  profiles.  During  oxide  growth  the  arsenic  is  rejected 
by  the  growing  oxide  region  [3],  and  instead  accumulates  inside  the  silicon 
rear  the  silicon-silicon  dioxide  interface.  The  accumulated  arsenic  can 
also  diffuse  from  the  interface  into  the  bulk  at  the  oxide  growth  temperature, 
and  the  two  mechanisms  of  segregation  and  diffusion  combine  to  produce  the 
final  arsenic  distribution.  At  900*"C  the  wot  oxidized  layer  grows  at  a 
faster  rate  than  the  arsenic  diffuses  away  from  the  Interface,  so  that  the 
entire  arsenic  dose  is  contained  within  a few  hundred  angstroms  of  the  surface. 
At  1100°C  the  situation  Is  reversed,  and  the  arsenic  diffuses  away  from  the 
interface  at  a faster  rate  than  the  silicon-silicon  dioxide  interface  advances. 
Since  the  oxide  layers  incorporate  the  entire  implanted  arsenic  distributions, 
it  is  no  surprise  that  the  final  arsenic  distributions  are  nearly  identical  for 
the  bare  surface  and  through-oxide  implanted  samples. 


15  -2  + 

Figure  4.13.  SIMS  results  for  5 x 10  cm  , 150  keV  As  implants 
following  2 hour  wet  oxidation  at  900°C.  Arsenic 
profile  (dots)  was  determined  after  the  Si02  layer 
was  chemically  removed;  on  the  bottom  the  arsenic 
profile  (squares)  was  determined  by  sputtering  through 
the  Si02  layer.  Apparent  in  the  lower  profiles  is  a 
decrease  in  depth  resolution  with  sputtering  time  and 
a loss  of  intensity  due  to  oxide  charging  effects. 
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In  examining  the  profiles  of  the  900°C  oxidized  samples,  we  notice 
the  limitation  of  the  SIMS  technique  for  near  surface  analysis.  Whereas  the 
arsenic  distribution  is  peaked  at  the  interface,  the  SIMS  technique  shows  a 
peak  located  away  from  the  interface,  due  to  a lack  of  chemical  equilibrium 
between  the  primary  ion  beam  and  the  substrate  for  the  initial  stages  of 
sputtering  [57].  This  effect  reduces  the  initial  secondary  ion  yield,  for 
the  reasons  discussed  in  section  4.2.3.  By  comparing  the  profiles  obtained 
from  sputtering  through  the  SiO^  layer  with  those  obtained  when  the  Si02  is 
chemically  removed  prior  to  examination  (Fig.  4.13),  we  again  see  profile 
broadening  and  reduction  in  ion  yield,  as  in  the  Auger  spectroscopy  results. 

Finally,  it  should  be  noted  that  one  anomaly  occurred  in  the  measure- 
ment of  the  sample  redistributed  by  thermal  oxidation  at  1100°C  following 
implantation.  The  profiles  obtained  for  samples  3 and  6 were  of  the  same 
shape,  but  the  absolute  intensity  for  sample  6 (bare  surface)  was  55%  larger 
than  that  obtained  from  sample  3 (through-oxide  implant).  It  has  previously 
been  observed  that  the  secondary  ion  yield  from  crystalline  samples  is  strongly 
orientation-dependent  [57],  and  it  appears  likely  that  for  sample  6 the  primary 
cesium  beam  in  the  SIMS  system  was  incident  upon  a preferred  crystalline 
direction.  However,  the  calibrate '-.i  for  sample  6 led  to  an  integrated  arsenic 
intensity  approximately  50%  greater  than  the  total  dose,  while  the  integrated 
intensity  for  sample  3 agrees  with  the  ion  dose  to  within  experimental 
uncertainty  (approximately  10%).  For  this  reason,  the  results  obtained  for 
sample  6 were  scaled  to  correspond  to  the  results  obtained  for  sample  3. 
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4.4.  Comparison  of  Techniques 

ihe  results  obtained  in  section  4.3  allow  us  to  make  several 
observations  about  the  relative  advantages  of  each  of  the  various  profiling 
techniques . 

Auger  spectroscopy  excels  in  the  measurement  of  near-surface 

features,  but  suffers  from  poor  impurity  detection  limits  (approximately 
20  -3 

4 X 10  cm  for  arsenic).  AES  has  excellent  lateral  resolution  [48]  (v;hich 

is  not  examined  in  this  study)  and  the  capability  to  monitor  more  than  one 

element  simultaneously.  The  GDOS  technique  is  well  suited  to  the  measurement 

of  implanted  profiles,  and  provides  easier  concentration  calibration  than 

any  of  the  other  techniques.  The  As  detection  limit  for  GDOS  is  between 

that  for  SIMS  and  AES.  By  using  multiple  analysis  and  detection  schemes  the 

GDOS  technique  can  simultaneously  monitor  the  presence  of  several  elements 

in  the  sample  simultaneously.  While  not  of  direct  concern  for  this  study, 

GDOS  is  also  the  least  expensive  depth  profiling  technique  of  the  three 

examined  in  this  study.  The  GDOS  technique  suffers,  however,  in  lack  of 

lateral  resolution  and  the  amount  of  sample  area  consumed  (approximately 

2 

0.5  cm  per  profile).  Also,  in  GDOS,  it  is  not  possible  to  examine  samples 
implanted  through  an  oxide  layer  without  first  chemically  removing  the  oxide 
layer.  Glow  discharge  optical  spectroscopy  is  a technique  still  in  its 
infancy,  and  as  more  work  is  put  into  the  development  of  this  technique  it 
should  be  possible  to  overcome  the  problems  encountered  in  profiling  the  annealed 
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Finally,  we  consider  SIMS.  lliis  technique  is  the  nost  sensitive 
and  also  the  most  expensive  of  the  techniques  employed  here.  The  SIMS 
system  used  in  this  study  detects  arsenic  concentrations  as  low  as  1 x cm 

using  a cesium  ion  source  [56],  although  a distinction  must  be  made  between 
ultimate  sensitivity  and  the  dynamic  range  of  the  instrument.  The  SIMS 
system  used  here  can  only  detect  changes  in  impurity  concentration  over 
three  orders  of  magnitude,  below  w'hich  impurities  sputtered  from  the  sides 
of  the  crater  walls  create  a background  ion  count  above  the  ultimate  detection 
limits  of  the  system.  The  SIMS  technique  excels  at  producing  the  shape  of 
the  dopant  vs.  depth  relationship;  however,  the  calibration  of  impurity 
concentrations  and  depth  scales  requires  care  and  considerable  effort. 
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5.  ELECTRICAL  MEASUREMENTS 

In  an  ion  implanted  sample  both  the  carrier  concentration  and 
carrier  mobility  vary  with  depth  from  the  sample  surface;  thus  any 
electrical  measurement  of  an  implanted  layer  represents  a weighted  average 
over  the  depth  distribution  of  implanted  ions.  By  combining  electrical 
measurements  with  successive  layer  removal,  it  is  possible  to  isolate  the 
contribution  from  each  layer  and  determine  its  average  carrier  concentration 
and  mobility.  By  comparing  such  measurements  on  samples  that  have  been 
implanted  through  thin  oxide  layers  with  samples  that  have  been  implanted 
into  bare  surfaces  under  the  same  conditions,  we  obtain  a direct  measurement 
of  the  effects  of  the  recoil  damage  on  the  electrical  properties  of  the 
implanted  layers.  In  addition,  these  measurements  allow  a direct  examination 
of  the  diffusion  and  activation  of  ion  implanted  arsenic  in  silicon. 

5.1.  Experimental  Procedure 

The  samples  used  for  this  part  of  the  study  w'ere  identical  to 
those  described  in  section  4.1  and  summarized  in  Table  111.  Prior  to  annealing, 
the  samples  were  cleaned  with  organic  solvents  and  H2O2 -based  cleaning 
solutions  [60].  All  anneals  performed  were  of  30  minutes  duration  under 
flowing  argon  gas.  Following  annealing,  the  samples  were  cleaved  into 
approximately  100  mil  squares,  cleaned,  and  then  etched  in  a buffered-HF 
solution.  Contacts  of  Au(l%Sb)  were  then  evaporated  onto  the  front  surface 
through  a shadow  mask  to  define  a van  der  Pauw  [61]  pattern.  These  contacts 
were  annealed  at  temperatures  between  400°C  and  425°C  (well  above  the  370° 

Au-Si  eutectic  [62])  for  approximately  20  seconds  and  then  quenched  in  distil  lec 
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water.  This  procedure  forms  alloy  junctions  that  penetrate  the  implanted  layers. 
Following  contact  formation,  the  samples  were  mounted  over  a different  shadow 
mask,  held  in  place  with  glycol  phthalate,  and  a groove  pattern  was  sandblasted 
into  the  surface  to  minimize  the  effects  of  finite  contact  size  [63].  In 
addition  to  providing  positive  location,  the  glycol  phthalate  also 
protected  the  sample  surface  during  the  sandb las t ing . The  completed  sample 
(Fig.  5.1)  was  next  removed  from  the  mask,  cleaned,  and  mounted  on  nylon  disks 
using  GK  7031  varnish  and  aligned  to  provide  proper  registration  in  the 
sample  holder  (Fig.  5.2).  F.lectricai  contact  to  the  sample  was  made  with 
spring  loaded  non-magnetic  Pogo  stick  contacts  made  by  the  Pylon  Co. 

Hall  effect  and  sheet  resistivity  measurements  are  made  using  a 
double  ac  Hall  effect  system,  described  elsewhere  [6u]  . in  this  method,  the 
magnet  is  excited  with  current  of  one  frequency  (250  Hz)  and  the  sample 
with  current  of  a different  frequency  (1  kHz).  Tne  Hall  voltage  then  appears 
at  the  difference  frequency.  Since  the  driving  frequencies  arc  phase  locked, 
one  measures  the  heterodyne  component  using  a lock-in  technique  as  shown 
in  the  schematic  diagram  (Fig.  5.3),  providing  excellent  noise  suppression. 

For  the  van  der  Pauw  style  double  ac  method,  the  sFieet  mobility  and  sheet 
carrier  concentration  are  given  by  [64] 

. n (’^ABCD  *’’^BCDA\  , /''ABCD')  _ ,, 
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Fig.  5.1.  Photograph  of  a typical  van  der  Pauw  sample  on  a nylon 
substrate.  The  contacts  are  isolated  by  sandblasted 
regions.  The  locating  slot  on  the  substrate  enables 
proper  positioning  of  the  sample  in  the  double  ac 
sample  holder. 
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Cross  sectional  view  of  the  sample  holder  which  is  placed 
inside  the  magnet  for  double  a -c  Hall  effect  measurements. 
The  inset  shows  an  enlarged  view  of  the  sample  and  spring 
mounted  contacts  in  the  holder  (after  McLevige  ^ £l_.  [64]) 
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If  the  contacts  are  labeled  clockwise,  represents  the 

voltage  developed  across  contacts  A and  B when  unit  current  is  passed 

through  contacts  C and  D.  The  other  symbols  are  defined  similarly. 

In  the  same  way,  AV  represents  the  voltage  measured  at  the  heterodyne 
BDAC 

frequency  across  contacts  B and  D for  unit  current  through  contacts  A and  C 

under  the  ac  magnetic  field  B(rms)  (202  gauss  for  the  present  system). 

R 

Here  f is  g tabulated  factor  given  by  van  der  Pauw  [61]  and  e is  the 

'*^BCDA 

electronic  charge.  By  successively  removing  layers  of  thickness  d^  the 

average  carrier  concentration  and  mobility  of  the  layer  can  be  determined 

by  measuring  the  sheet  resistivity  o and  sheet  Hall  coefficient  R before 

s.  s. 

stripping  and  the  same  quantities  after  stripping  (i.e.,  D ,R  )•  The 

Sj-1 

average  mobility  in  the  layer  is  given  by 
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and  the  average  carrier  concentration  is  given  by  [27] 


1/0  -I/03 

= j-1 

"j  ed.u. 
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The  derivation  of  these  equations  based  on  an  equivalent  circuit 
model  is  presented  in  the  Appendix;  however,  these  can  also  be  derived  from 
the  Boltzmann  transport  equation  [65].  In  the  derivation  of  these  equations 
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In  tilt  same  wa\  , represents  the  voltage  measured  at  the  ht  cerodyne 

frequency  across  contacts  l>  and  D iiir  unit  turrent  througli  contacts  A and  C 
under  the  tic  magnetic  field  li(rms)  (202  gauss  for  tlie  present  system). 

[U  Lc  it;  a tabulattd  factor  given  by  van  der  Fanw  |61J  and  e is  the 
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electronic  charge.  By  successively  removing  laycr.s  of  thickness  d.,  the 
average  carrier  concentration  and  mobility  of  the  layer  can  be  determined 
by  measuring  the  sheet  re.sistivity  o and  sheet  Hall  coefficient  R before 

Sj  Sj 

stripping  and  the  same  quantities  after  stripping  (i.e.,  D,  ,R  ) The 
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average  mobility  in  the  layer  is  given  by 
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The  derivation  of  these  equations  based  on  an  equivalent  circuit 
model  is  pre.scnted  in  the  Appendix;  however,  liiese  can  also  be  derived  from 
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it  is  assumed  that  the  conductivity  mobility  and  the  Hall  mobility  arc  equal. 

Ihe  two  mobilities  are  equal  in  the  case  of  impurity  banding  [66],  which 

18  - 3 

takes  place  in  silicon  for  donor  concentrations  of  3 x 10  cm  [67].  Thus, 

18  “3 

for  concentrations  greater  than  3 x 10  cm  , the  carrier  concentrations 
will  be  correct  to  within  the  experimental  limitations  of  the  technique. 

Following  the  measurement  of  surface  sheet  resistivity  and  sheet 
Hall  coefficient,  the  successive  layer  removal  is  accomplished  by  chemical 
etching.  Before  the  etch  is  begun,  Apiezon  W (black  wax)  docs  of  approximately 
0.5  mm  are  melted  onto  the  sample  surface  outside  the  corners  of  the  contact 
pads.  These  dots  serve  Co  protect  the  original  surface  and  allow  for  a direct 
depth  calibration  through  the  use  of  a surface  prof i lometer . 

Controlled  etching  is  performed  by  immersing  the  sample  in  a 
beaker  partly  filled  with  the  etch  solution  for  a fixed  period  of  time.  The 
mounted  sample  is  held  with  Teflon-coated  tweezers  and  the  solution  is  manually 
stirred  at  approximately  120  rpm  during  the  etching  process.  Following  the 
etching,  the  sample  is  immediately  placed  in  water  and  then  dried  with 
compressed  dry  nitrogen.  The  sheet  resistivity  and  Hall  measurements  are 
taken,  and  Che  entire  process  is  repeated  suc.cessive ly  until  the  implanted 
layer  is  removed. 

The  standard  silicon  planar  etch  [68]  is  300  ml  nitric  acid,  100  ml 
acetic  acid,  4 ml  hydrofluoric  acid;  however,  its  etch  rate  was  determined 
experimentally  to  be  approximately  5 wn/minute.  Since  Che  concentration  of 
hydrofluoric  acid  controls  the  etch  rate  [69],  a more  dilute  etch  was 
constructed  by  mixing  10.4  ml  of  a mixture  consisting  of  100  ml  acetic  acid 
and  4 ml  of  hydrofluoric  acid  with  90  ml  of  acetic  acid  and  300  ml  of  nitric 
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acid.  The  more  dilute  etch  was  used  during  this  studv.  Itie  etch  rate  varied 
from  one  mixture  to  another  from  a minimum  of  230  X/minu'.e  to  a maximum  of 
580  X/minute,  but  was  fairly  consistent  within  a give-  batch.  Tlie  variation 
in  etch  rate  from  one  solution  to  the  other  indicate?  the  difficult-.-  in 
controlling  such  small  concentrations  of  n ,-drof  luoric  a.-.id  in  the  solutions. 

It  is  also  possible  that  other  factors,  .s-jch  as  extent  of  crystalline  damage 
and  the  position  of  the  Fermi  level  in  the  sample,  may  also  have  had  an 
effect  on  the  etch  rate. 

The  concentration  dependence  of  tr.e  etch  was  deterrriined  by  etching 

through  a sample  with  a nearly  gaussian  impurity  distribution,  measured  by 

the  differential  Hall  effect  method.  Tliis  sample  was  etched  past  the  peak 
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of  the  imparity  distribution  until  the  impurity  cuncentrat-Ion  (£.3  x 10  cm  ) 
was  7.57,  of  the  peak  concentration.  Tne  dots  were  then  removed,  other  dots 
formed  in  different  areas  of  the  sample,  and  an  equal  number  of  steps  remov'ed . 
The  step  heights  were  then  measured  with  a Sloan  Dektak.  The  etch  race 
thrc'ugh  the  region  with  the  maximum  imp-„rity  cur.centratic:.  was  300  + 30  X/mirute 
while  the  etch  rate  for  the  lower  concentration  regicr"  was  290  + 29  ^/minute. 
Measurements  on  similar  samples  using  a different  batch  of  etch  mixture 
yielded  332  + 30  X/minute  for  the  peak  of  the  distribution  and  325  + 30  X/minute 
for  the  tail  region  of  the  distribution.  For  all  the  samples  examined  the 
change  in  etch  rate  with  concentration  was  less  than  the  uncertainty  in  the 
etch  rate.  The  samples  which  were  steam  oxidized  following  implantation  showed 
the  greatest  variations.  For  such  samples,  the  etch  rate  was  typically 
380  + 40  ^/minute  for  the  peak  of  the  distribution  and  320  + 30  S/minute  for 
the  tail  regions  of  the  distribution.  In  all  cases,  the  thickness  of  the  layer 
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removed  was  calculated  from  the  measured  etcli  rate  for  that  number  of  steps 
and  not  from  an  etch  rate  averaged  over  the  total  profile. 

5.2.  Sheet  Measurements 

The  formation  of  emitters  for  silicon  bipolar  transistors  is  one  of 

the  major  applications  of  arsenic  implantation  in  silicon.  Since  the  emitter 

sheet  resistivity  plays  an  important  part  in  the  determination  of  transistor 

gain  [70] , measurement  of  the  sheet  resistivity  of  arsenic  implanted  layers 

has  immediate  technical  importance.  In  addition,  these  measurements,  when 

combined  with  measurements  of  sheet  Hall  coefficient,  yield  a direct  estimate 

of  the  activation  of  the  implanted  ions  for  the  bare  surface  implants. 

Difficulties  are  encountered,  however,  in  estimating  the  activation  of  the 

through-oxide  implants  since  not  all  of  the  total  dose  reaches  the  silicon 

substrate.  For  those  cases,  one  must  compare  the  integrated  electrically 

active  concentration  to  the  integrated  total  impurity  concentration  to  obtain 

a true  value  for  activation.  In  this  study,  all  values  quoted  represent  ar. 

average  of  measurements  from  at  least  four  identical  samples. 

Sheet  resistivity  measurements  have  recently  been  performed  through- 

oxide  implanted  samples  [71],  but  these  were  only  for  1000°C  anneals. 

No  mention  was  made  of  electrical  activation  and  no  measurements  of  carrier 

or  mobility  profiles  were  attempted.  The  values  obtained  in 

^ . udy  for  the  .sheet  resistivities  and  sheet  Hall  c oe  f f i c i onl  !■. 

of  the  bare  surface  and  through-oxide  implants  following  600°C  annealing 

are  summarized  in  Table  IV.  From  these  results  it  can  be  seen  that  the 
15  -2 

1 X 10  cm  implants  attain  almost  identical  sheet  resistivities  following 
600°  annealing.  It  can  also  be  seen  that  the  1 x lO^^cm  ^ bare  surface  implant 
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TABLE  IV 


Average  sheet  resistivities  and  average  sheet  carrier  concentrations  for 
As^  implanted  ( 100)  p-type  silicon  substrates  following  600°C  30  minute 
annea ling . 


Sample 

Type 

As 

Implant 

Surface  condition 
before  Implantation 

Average 

Sheet 

Resistivity 

Average  Sheet 
Carrier 
Concentration 

1 

15  -2 

5x10  cm  150  keV 

150  X SiO^ 

65  n/D 

2 . IxlO^^cm'^ 

5 

5xl0^^cm’^150  keV 

Bare 

56  !!/□ 

15  -2 

1.9x10  cm 

13 

1x10 '"^cm  ^200  keV 

840 X 3 0^ 

98  Q/G 

g „ 14  -2 

8.2x10  cm 

14 

lxlO^^cm’^200  keV 

Bare 

96  Q/n 

14  -2 

9.9x10  cm 

h;is  readied  nearly  complete  electrical  activity  following  the  600  C ameal<-i 


consistent  w’ith  the  work  of  Crowder 


activity  for  peak  impurity  concentrations  less  than  2 x iU  cn. 


snows  approxima 


ions,  which  is  reasonable  in  vie 


impurity  concentration  of  7 x 10“  cm 


activity,  the>  were  not  subjected  to  further  anneals.  Samples  subjected  ti 


s were  examinea 


‘lOO'C  or  1000  C,  and  the  results  are  summarized  in  Fig.  5.6.  Again,  i 


emperaturc  to  1 


the  through  oxide  implanted  samples  compared  to  the  bare  surface  impla 


in  view  of  the  vast  crystallographic  damage  seen  by  trar.smi.s? ion  electron 


microscopy  [4,7]  and  by  He  ion  backscattering  [5]  for  through  oxide  implanted 


samples  annealed  at  1000  C.  When  the  5 x 10  cm  implanted  samples  are  a 


for  5 minutes  in  dry  nitrogen,  followed  by  SiO.  thermal  growth  which  incorporates 


the  entire  implanted  layer,  sheet  resistivities  and  carrier  concentrations  were 


nearly  identical  for  the  bare  surface  and  through  oxide  implants  (Table  V) 


Figure  5.4.  Variation  in  sheet  carrier  concentrations  (top  graph) 

and  sheet  resistivities  (bottom)  with  annealing 

15-2 

temperature  for  5 x 10  cm  150  keV  As  implants  into 
bare  surfaces  (squares)  and  through  approximately  150 
Si02  (circles). 


84 


TABL8  V 

Average  sheet  resistivities  and  average  sheet  carrier  ccncenrrat tons  for 
15-2  + 

5 X 10  cm  150  keV  As  -implanted  ^ lOO)  p-type  silicon  substrates 
following  redistribution  by  thermal  oxidation.  Oxide  tbickr.ess  was  deter- 
mined by  e 1 li psomet ry . 


Samp le 
Type 

Surface  before 
Implantation 

Oxidation 
Tempera  ture 

Final  Si02 
Thickness 

Sheet 

Resistivity 

Ca  rr  le  r 
Concent  ratio 

2 

~150S  SiO^ 

900°C 

32  72  R 

126  I/O 

- 0 - 
/ . 3x lu  cm 

3 

~ 150S  SiO^ 

1100°C 

4432  S 

20  0/D 

4 . / xlO  ^^crr. 

4 

Bare 

900°C 

3338  X 

128  Q/n 

7 . lx  10  ^ "'em 

6 

Bare 

1100°C 

4382  S 

20  -UU 

4 .6xl0^'Cn 
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also  is  surprising,  in  view  of  the  fact  that  the  2 to  1 increase  in  volume  in 
the  surface  silicon  layer  due  to  thermal  oxidation  is  known  to  propagate 
surface  defects  further  into  the  silicon  substrate  [73]. 

Finally,  an  attempt  was  made  to  determine  an  activation  energy 
for  the  annealing  of  sheet  resistivity  and  sheet  carrier  concentration  of  the 
5 X lO^^cm  ^ implants.  Isothermal  anneals  were  performed  at  600°C  for 
30  minutes,  60  minutes,  and  150  minutes,  but  no  change  in  sheet  carrier  concen- 
tration nor  sheet  resistivity  could  be  observed  among  the  various  anneals, 
within  experimental  uncertainty.  Apparently,  annealing  operative  at  600°C 
occurs  rapidly  (<  30  rain). 


5.3.  Carrier  Concentration  and  Mobility  Profiles 

The  only  way  to  obtain  a true  estimate  of  the  impurity  activation 
and  carrier  transport  properties  of  the  through-oxide  implanted  samples  is  by 
combining  differential  Hall  effect  measurements  with  a technique  that  measures 
the  total  impurity  distribution,  such  as  SIMS.  These  mea suremer, ts  are 
performed  on  the  through-oxide  implanted  samples  and  also  on  the  equivalent 
bare  surface  implanted  wafers.  This  comparison  provides  a direct  measurement 
of  the  effects  of  the  damage  produced  by  recoil  implanted  oxygen  on  the 
transport  properties  of  through-oxide  implanted  layers.  Although  the  inter- 
pretation of  the  results  obtained  is  complicated  by  the  fact  that  the  grown 
oxides  were  of  nonuniform  thickness,  we  shall  see  that  the  depth  profiles 
verify  the  trends  seen  from  the  sheet  resistivity  measurements  and  provide  a 
partial  explanation  of  these  trends. 

In  Fig.  5.5  we  present  the  carrier  concentration  and  mobility  profile 
15  -2 

for  the  5 X 10  cm  150  keV  As  implant.  Comparing  the  electrical  profile 
with  the  total  imnurity  profile  indicates  a concentration-dependent  activation 
after  600°C  annealing  consistent  with  the  measurements  of  Crowder  [72]. 
Integrating  the  carrier  concentration  over  the  depth  distribution  yields  a 
result  in  excellent  agreement  with  the  sheet  carrier  concentration  reported 
in  Table  IV.  In  the  plot  of  the  mobility  as  a function  of  depth,  the 
experimentally-obtained  mobilities  are  compared  to  mobility  values  expected 
for  the  experimentally-determined  carrier  concentrations.  The  expected 
values  of  mobility  for  the  given  carrier  concentrations  are  obtained  from  an 
eleven-term  polynomial  fit  to  the  empirical  resistivity-carrier  concentration 
relationship  for  arsenic-doped  silicon  determined  by  Fair  and  Tsai  [59]. 


Figure  5.5.  Carrier  concentration  and  mobility  profiles  for  600°C 

15  -2  + 

30  minute  anneal  of  5 x 10  cm  150  keV  As  implants 
into  bare  Si  surface.  Top:  comparison  of  electrically 

active  arsenic  concentration  obtained  by  differential 
Hall  effect  measurements  (circles)  with  total  arsenic 
concentration  obtained  by  SIMS  (solid  line).  Bottom: 
comparison  of  mobility  values  obtained  by  differential 
Hall  effect  measurements  with  values  expected  from 
polynomial  fit  of  Fair  and  Tsai  [59]  to  empirical 
carrier  concentration-mobility  relationship  for  As-doped 
Si. 


Carrier  Mobility  icm/v-secj  Cone 


It  is  seen  that  the  experimentally  observed  mobilities  agree  well  with  the 
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expected  values.  Annealing  the  5 x 10  cm  150  keV  As  implant  to  1000  C 
(Fig.  5.6)  produces  excellent  electrical  activation,  in  good  agreement  with 
the  other  studies  of  implanted  arsenic  [59],  and  the  experimentally-determined 
mobilities  exhibit  excellent  agreement  with  the  expected  mobilities. 
Integrating  the  electrically  active  profile  yields  results  in  excellent  agree- 
ment with  the  sheet  carrier  concentration  values.  Throughout  this  study, 
the  integrated  carrier  concentrations  never  differed  from  the  measured  sheet 
carrier  concentrations  by  more  than  10%,  and  thus  the  agreement  will  no 
longer  be  noted  in  the  descriptions  of  the  remaining  profiles. 

In  Figs.  5.7  and  5.8,  we  examine  the  corresponding  samples  implanted 
through  thin  oxides.  Since  the  oxide  thickness  for  the  implanted  layer  was 
found  to  be  nonuniform,  difficulties  are  encountered  in  correlating  the  depth 
scales  between  the  electrically-active  profile  and  the  SIMS  profile. 

Throughout  this  study,  effort  was  made  to  take  samples  for  elemental 
profiling  and  samples  used  for  electrical  measurements  from  adjacent  areas 
on  the  implanted  wafer.  Despite  these  efforts,  examination  of  Fig.  5.7 
indicates  that  the  oxide  thickness  for  the  area  of  the  wafer  used  for  SIMS 
profiling  was  slightly  thicker  than  the  oxide  thickness  for  the  area  of  the 
wafer  used  for  electrical  profiling.  Thus  the  depth  distributions  peak  at 
different  values  from  the  silicon-silicon  dioxide  interface.  Nevertheless,  we 
observe  that  electrical  activation  similar  to  the  bare  surface  implant  is 
attained,  and  that  surprisingly,  the  mobilities  obtained  for  the  through-oxide 
implants  agree  well  with  the  expected  values  for  arsenic-doped  silicon. 
Annealing  the  through  oxide  implanted  sample  to  1000°C  - where  the  residual 


Figure  5.6.  Carrier  concentration  and  mobility  profiles  for  1000°C 
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30  minute  anneal  of  5 x 10  cm  150  keV  As  implants  into 
bare  Si  surface.  Top:  comparison  of  electrically  active 

arsenic  concentration  obtained  by  differential  Hall  effect 
measurements  (circles)  with  total  arsenic  concentration 
obtained  by  SIMS  (solid  line) . Bottom:  comparison  of 

mobility  values  obtained  by  differential  Hall  effect 
measurements  with  values  expected  from  polynomial  fit  of 
Fair  and  Tsai  [59]  to  empirical  carrier  concentration- 
mobility  relationship  for  As-doped  Si. 


Figure  5.7.  Carrier  concentration  and  mobility  profiles  for  600°C 

15  -2  + 

30  minute  anneal  of  5 x 10  cm  150  keV  As 
Implants  through  approximately  150  X Si02 • 

Top:  comparison  of  electrically  active  arsenic 

concentration  obtained  by  differential  Hall  effect 
measurements  (circles)  with  total  arsenic  concentra- 
tion obtained  by  SIMS  (solid  line).  Bottom: 
comparison  of  mobility  values  obtained  by  differential 
Hall  effect  measurements  with  values  expected  from 
polynomial  fit  of  Fair  and  Tsai  [59]  to  empirical 
carrier  concentration-mobility  relationship  for 
As -doped  Si. 


o 

Figure  5.8.  Carrier  concentration  and  mobility  profiles  for  1000  C 

15  -2  + 

30  minute  anneal  of  5 x 10  cm  150  keV  As  implants 
through  approximately  150  X Si02 • Top:  comparison 

of  electrically  active  arsenic  concentration  obtained 
by  differential  Hall  effect  measurements  (circles) 
with  total  arsenic  concentration  obtained  by  SIMS 
(solid  line).  Bottom:  comparison  of  mobility  values 

obtained  by  differential  Hall  effect  measurements  with 
values  expected  from  polynomial  fit  of  Fair  and  Tsai  [59J 
to  empirical  carrier  concentration-mobility  relationship 
for  As-doped  Si . 
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' damage  was  originally  observed  [4]  - again  shows  electrical  activation  identi- 

cal to  the  bare  surface  implant  and  surprisingly  good  mobility  values.  Even 
more  surprisingly,  when  silicon  dioxide  layers  arc  thermally  grown  to  incorporate 
the  entire  implanted  layer,  the  carrier  and  mobility  profiles  for  the  bare 
^ surface  and  through  oxide  implants  are  identical,  to  within  experimental 

uncertainties.  Ihus  in  Figs.  5.9  and  5.10,  only  one  set  of  points  is  shown 
to  indicate  the  profiles  obtained  from  both  implant  conditions.  In  Fig.  5.9, 
we  again  see  that  for  the  1100°C  redistribution  nearly  complete  electrical 
activation  is  obtained,  and  that  experimental  mobilities  agree  well  with  the 
expected  values.  Apparently  the  diffusion  coefficient  of  arsenic  in  silicon 
at  1100°C  is  large  enough  that  the  arsenic  displaced  by  the  moving  silicon- 
silicon  dioxide  boundary  diffuses  away  from  the  interface  fast  enough  to 
prevent  an  appreciable  impurity  pile-up.  Such  is  not  the  case,  however,  for 
the  900°C  oxidation.  Here  the  arsenic  has  accumulated  within  a few  hundred 
angstroms  of  the  silicon-silicon  dioxide  interface.  Such  shallow  profiles  are 
extremely  difficult  to  measure  accurately.  The  chemical  etch  rate  was 
appreciably  faster  for  the  peak  of  this  distribution  than  for  the  tail  region 
(section  5.1)  while  the  SIMS  data  overestimates  the  depth  of  the  distribution 

^ due  to  matrix  effects  and  differential  sputtering.  These  problems  make  it 

difficult  to  correlate  the  e lectrica lly-active  concentration  profile  with  the 

I 

elemental  profile,  although  both  methods  gave  consistent  readings  within  them- 

( 

selves.  Good  mobility  and  activation  are  also  obtained  for  annealing  of 
samples  implanted  through" thicker  oxides  (Fig.  5.11  - 5.12). 

The  work  reported  in  this  chapter  represents  the  first  complete  set  " 
of  electrical  measurements  of  the  transport  properties  of  through -oxide 

I 

I 

1 

J 


1 


•“'I  — 
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Figure  5.  9.  Carrier  concentration  and  mobility  profiles  following 
redistribution  by  wet  oxidation  for  60  minutes  at 
llOO^C  of  5 X lO^^cm  ^ 150  keV  As^  implants.  Results 
for  bare  surface  and  through  oxide  implants  identical. 
Top:  comparison  of  electrically  active  arsenic 

concentration  obtained  by  differential  Hall  effect 
measurements  (circles)  with  total  arsenic  concentration 
obtained  by  SIMS  (solid  line).  Bottom:  comparison  of 

mobility  values  obtained  by  differential  Hall  effect 
measurements  with  values  expected  from  polynomial  fit 
of  Fair  and  Tsai  [59]  to  empirical  carrier  concentration- 
mobility  relationship  for  As-doped  Si. 


'T 
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Figure  5.10.  Carrier  concentration  and  mobility  profiles  following 

redistribution  by  wet  oxidation  for  120  minutes  at  900°C 
15  -2  + 

of  5 X 10  cm  150  keV  As  implants.  Results  for  bare 
surface  and  through  oxide  implants  identical.  Top: 
comparison  of  electrically  active  arsenic  concentration 
obtained  by  differential  Hall  effect  measurements  (circles) 
with  total  arsenic  concentration  obtained  by  SIMS  (solid  line). 


Figure  5.11.  Carrier  concentration  and  mobility  profiles  for  600  C 

15  -2  + 

30  minute  anneal  of  1 x 10  cm  200  keV  As  implants 
through  approximately  840  X SIO2.  Top:  comparison 

of  electrically  active  arsenic  concentration  obtained 
by  differential  Hall  effect  measurements  (circles) 
with  total  arsenic  concentration  obtained  by  SIMS 
(solid  line). 


Figure  5.12.  Carrier  concentration  and  mobility  profiles  for  1000°c 

15  -2  + 

30  minute  anneal  of  1 x 10  cm  200  keV  As  implants 
through  approximately  840  X Si02*  Top:  comparison 

of  electrically  active  arsenic  concentration  obtained  by 
differential  Hall  effect  measurements  (circles)  with 
total  arsenic  concentration  obtained  by  SIMS  (solid  line) 
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implanted  silicon.  It  is  seen  that  good  electrical  properties  of  the  through- 
oxide  implanted  samples  suggested  by  the  sheet  resistivity  measurements  are 
due  to  two  factors;  activation  for  the  through-oxide  implant  identical  to 
bare  surface  implanted  silicon  and  the  attainment  of  carrier  mobilities  in 
the  through  oxide  implanted  layers  identical  to  mobilities  obtained  for 
bulk  silicon.  These  results  could  hardly  have  been  suspected  prior  to  this 
investigation,  due  to  the  known  complexing  of  arsenic  and  oxygen  in  radiation- 
damaged  silicon  [8]  and  the  extensive  damage  see  for  through -oxide  implanted 
samples  [4,5,7].  These  results  do,  however,  provide  for  a possible  inter- 
pretation, which  we  shall  develop  in  the  next  chapter. 

5.4.  Discussion  of  Uncertainties 

The  major  sources  of  uncertainty  in  determining  the  electrically 
active  impurity  profile  and  the  mobility  profile  for  the  implanted  layer  are 
nonuniformity  of  layer  thickness,  uncertainty  in  measurement  of  sample  currents 
and  voltages,  uncertainty  in  magnetic  field  strength,  and  junction  leakage 
fr-m  the  implanted  layer.  From  Eq.  (5.1),  it  is  seen  that  the  uncertainties 
in  sheet  Hall  coefficient  and  sheet  resistivity  are  directly  proportional 
to  the  uncertainty  in  the  measurement  of  sample  voltages,  while  the  uncertainty 
in  the  sheet  Hall  coefficient  is  inversely  proportional  to  both  the  magnetic 
field  strength  and  the  sample  current.  For  the  present  system,  the  magnetic 
field  is  202  + 2 gauss,  which  would  create  a 17o  uncertainty  in  the  sheet  Hall 
coefficient.  The  sample  current  is  determined  by  measuring  the  voltage 
developed  across  a 100  precision  resistor  using  a PAR  HR-8  lock-in  amplifier, 
which  is  read  out  on  a Keithley  160  digital  multimeter.  Sample  voltages  are 
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measured  directly  using  a PAR  186A  lock-in  amplifier  and  are  read  out  on  a 
Digitec  dc  digital  voltmeter.  Drift  of  the  reference  frequency  for  the 
HR-8,  voltage  loss  in  the  connecting  cables,  variations  in  contact  pressure 
for  the  Pogo  stick  contacts,  and  electrical  noise  all  combine  to  produce 
less  than  a 3%  spread  in  the  measured  voltage?  foi  a given  sample. 

To  understand  how  the  uncertainties  in  the  sheet  resistivity  and 
sheet  Hall  coefficient  affect  the  experimentally  determined  carrier 
concentrations  and  mobilities,  one  must  consider  the  mathematical  form  of 
Eqs . (5.4)  and  (5.5).  This  can  best  illustrated  by  first  considering  the 
effects  of  layer  removal  on  a uniformly  doped  sample,  and  then  extending  this 
analysis  to  the  experimentally  encountered  profiles. 

For  a uniformly  doped  sample,  graphs  of  sheet  resistivity  and  sheet 

Hall  coefficient  with  depth  form  hyperbolas,  the  two  quantities  being  related 

by  a constant.  The  above  can  be  verified  by  examination  of  Eqs.  (A.  17)  and 

(A.  18)  of  the  Appendix  for  n^  and  constant.  If  we  now  consider  an  actual 

impurity  profile  with  a slight  slope  (Fig.  5.13),  we  see  that  the  sheet 

resistivity  and  sheet  Hall  coefficients  are  no  longer  related  by  a constant, 

but  rather  that  the  sheet  Hall  coefficient  increases  more  rapidly  with 

depth  than  the  sheet  resistivity  (Fig.  5.14).  This  behavior  is  expected 

since  the  sheet  Hall  coefficient  is  essentially  proportional  to 

while  the  sheet  resistivity  is  proportional  to  1/n  where  n and  u 

^ eff  eff  etf  etf 

are  "average"  carrier  concentrations  and  mobilities,  respectively.  Because 
lower  carrier  concentrations  have  higher  mobilities,  the  loss  of  free  carriers 
in  the  sheet  resistivity  is  partly  offset  by  their  increased  mobility,  while 
this  effect  does  not  occur  for  the  sheet  Hall  coefficient.  The  changes  in 


Figure  5.13.  Measured  depth  variation  of  carrier  concentration 

and  mobility  for  nearly  uniform  profile  with  slight 
gradient  (1000°C  30  minute  anneal  of  the  5 x 10 
150  keV  As  implant  through  approximately  150  X of  Si02)- 
Scatter  in  the  first  two  points  is  due  to  the  small 
percentage  change  in  sheet  resistivity  at  the  surface 
and  the  scatter  near  1000°C  is  due  to  a bad  data  point 
in  the  sheet  measurements  at  1000 


Sheet  Hall  Coefficient  (cm  /coul) 
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Figure  5.14.  Measured  variation  in  sheet  resistivity  and  sheet  Hall 

coefficient  with  depth  for  the  sample  shown  in  Fig.  5.13 

15  "2 

(1000°C  30  minute  anneal  of  5 x 10  cm  150  keV  As 
implant  through  Si02)* 


Sheet  Resistivity  (ohms/square) 
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sheet  resistivity  and  sheet  Hall  coefficients  are  oven  more  drastic  for  the 
case  of  a gaussian  impurity  distribution  (J'igs.  5.J5  and  5.1b). 

ihe  graphs  of  tlie  raw  data  presented  in  i'rgs.  5.13  and  5.14  help 
to  explain  otlier  limitations  of  the  techr.ir  e . iie.a  r the  snrlace,  i.he  sheet 
resistivity  suffers  only  a small  percentage  change  on  the  removal  of  the  first 
few  layers  (Fig.  5.13),  and  thus  uncertainties  in  the  measurement  of  each 
sheet  resistivity  can  amount  to  a significant  percentage  of  the  difference 
between  successive  sheet  resistivities.  Because  Eq . (5.4)  requires  the 
difference  between  the  inverse  sheet  resistivities  and  also  the  difference  of 
their  squares,  small  uncertainties  in  the  measurement  of  the  sheet 
resistivity  (e.g.  37,)  can  yield  large  variations  (approximately  507,)  in  sheet 
mobility  and  carrier  concentration  near  the  surface,  as  is  seen  in  the  profile 
for  the  1000°C  anneal  (Fig.  5.14).  Since  Eqs . (5.4)  and  (5.5)  rely  on  the 
difference  between  two  sets  of  data  points,  one  bad  measurement  affects  two 
sets  of  data  points:  the  values  obtained  by  comparing  the  bad  data  point 

with  the  preceding  step  and  the  values  obtained  by  comparing  the  bad  data 
point  with  the  following  step.  Thus,  a bad  data  point  at  1000  for  the 
1000°C  annealed  sample  shown  in  Fig.  5.13  leads  to  incorrect  values  of 
mobility  and  carrier  concentration  at  800  ^ and  1200  R (Fig.  5.14).  Even  at 
the  peak  of  the  gaussian  distribution  where  the  change  in  both  sheet  resistivity 
and  sheet  Hall  coefficient  is  greatest  and  thus  the  percentage  uncertainty 
in  the  change  is  the  least,  values  of  carrier  concentration  and  mobility  are 
extremely  sensitive  to  variations  in  sheet  resistivity,  as  can  be  seen  from 
Table  VI.  For  this  reason,  the  effect  of  junction  leakage  on  these  measure- 
ments is  of  extreme  concern.  As  a rule  of  thumb,  the  data  is  no  longer 


Ill 


Figure  5.15.  Measured  depth  variation  of  carrier  concentration  and 
mobility  for  sample  with  gaussian-type  impurity 
distribution  (600°C  30  minute  anneal  of  5 x lO^^cm 
150  keV  As^  implant  through  approximately  150  X Si02). 
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Figure  5.16.  Measured  variation  in  sheet  resistivity  and  sheet  Hall 

coefficient  with  depth  for  sample  shown  in  Fig.  5.15 

15  -2 

(600°C  30  minute  anneal  of  5 x 10  cm  150  keV  As 
implant  through  SiO-). 
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considered  valid  when  the  maximum  sample  current  through  the  implanted  layer 
is  less  than  fifty  times  the  leakage  current. 

While  the  differential  Hall  effect  method  is  extremely  sensitive  to 
the  value  of  sheet  resistivity,  it  is  not  strongly  affected  by  the  thickness 
of  the  layer  removed.  Layer  thickness  does  not  enter  into  the  calculation 
of  average  mobility  of  Eq . (5.4),  and  enters  inversely  into  the  c.Tlculation 
of  carrier  concentration  from  Eq . (5.5).  On  the  other  hand,  layer  thickness 
enters  directly  into  the  depth  calibration.  In  this  study,  etching  was 
performed  for  30+2  seconds/step,  giving  a minimum  7%  variation  in  the 
uniformity  of  the  etch  step,  disregarding  other  factors  possibly  affecting 
the  etch  rate,  such  as  residual  damage  or  carrier  concentration  dependence 
of  the  etch.  Thus,  in  the  construction  of  the  profiles  presented  in  this 
chapter,  the  7%  uncertainty  in  carrier  concentration  due  to  the  uncertainty 
in  layer  uniformity  was  accepted  rather  than  trying  to  estimate  "smoothed" 
values  of  resistivity  or  Hall  coefficient;  errors  in  smoothing  these  values 
would  lead  to  disproportionately  larger  errors  in  carrier  concentration  or 
mobility.  The  carrier  concentration  and  mobility  profiles  represent  the 
overlay  of  from  three  to  five  measurements  on  identical  samples  with  the  range 
of  values  obtained  listed  as  the  uncertainty  in  the  measurements. 


j 
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6.  DISCUSSION  AND  CONCLUSIONS 


It  is  well  understood  that  when  bare  silicon  surfaces  are  implanted 

14  -2 

with  arsenic  to  doses  greater  than  10  cm  , an  amorphous  region  is  formed 

at  the  sample  surface  due  to  the  extensive  damage  created  by  the  energetic 

arsenic  ions  [27].  For  annealing  temperatures  greater  than  550°C,  this 

amorphous  layer  regrows  epitaxially  on  the  undamaged  substrate  [72].  For 
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peak  arsenic  concentrations  less  than  approximately  2.0  x 10  cm  , all 
the  implanted  arsenic  ions  become  electrically  active  [72],  as  is  seen  in 
this  study.  Implantations  with  higher  peak  arsenic  concentrations  result 
in  incomplete  electrical  activation,  the  extent  of  which  depends  on  the 
annealing  time  and  temperature  used  [59] . This  incomplete  activation  is 
mainly  due  to  the  formation  of  electrically-inactive  arsenic  clusters  for 
high  dopant  concentrations,  also  observed  in  silicon  grown  with  high 
arsenic  content  [74].  The  percent  activation  is  improved  with  higher  annealing 
temperatures,  due  to  the  lowering  of  the  peak  arsenic  concentration  due  to 
diffusion.  Support  for  this  interpretation  is  provided  by  the  electrical 
results  for  the  bare  surface  implants  reported  here  (Chapter  3'.  Tie 
600  C anneal  of  the  5 x lO^^cm  ^ bare  surface  (Fig.  5.5)  results  in  the 
lowest  percentage  activation  at  the  peak  of  the  arsenic  concentration,  while 
the  regions  of  the  profile  with  impurity  concentrations  below  a pproximate ly 
1.5  X 10^*^cm  ^ show  complete  electrical  activation  within  the  experimental 
uncertainty.  As  the  annealing  temperature  approaches  1000°C,  the  profile 


has  spread  such  that  most  of  the  impurity  concentration  is  less  than 
20  -3 

2 X 10  cm  (Fig.  5.6)  and  the  sample  shows  90  to  95%  electrical  activation. 
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Consistent  with  this  interpretation,  the  activation  of  the  5 x lO^^cm'^ 

bare  surface  implanted  sample  after  .recrystallization  is  constant  for  annealing 

times  up  to  150  minutes  (Section  5.2).  This  correlates  well  with  the  fact 

that  measurable  arsenic  diffusion  has  not  been  observed  for  600°C  annealing  [1]. 

14  -2 

From  these  results  a commonly  accepted  view  of  high  dose  (>  10  cm  ) 

As  implantations  into  bare  surfaces  has  been  proposed.  According  to  this 
model,  the  implantation  creates  an  amorphous  layer  on  top  of  an  essentially 
undamaged  substrate.  This  layer  regrows  during  annealing  resulting  in 
arsenic  activation  similar  to  that  observed  in  arsenic-doped  silicon  boules. 

Thus  it  is  suggested  that  the  eptiaxial  regrowth  of  bare  surface  implanted 
amorphous  samples  is  a simple  example  of  an  amorphous -to -crys ta 1 line  phase 
transformation  [75].  Channeling  effect  measurements  on  bare  surface  implanted 
samples  indicate  a linear  regrowth  rate  and  a well  defined  activation  energy 
[75] . Both  the  quality  of  the  regrown  layer  and  the  regrowth  rate  are 
functions  of  substrate  orientation  [6,76],  with  the  most  disorder  seen  for 
( 11 l) -oriented  substrates  [6]  with  damage  consisting  of  extensive  dislocation 
networks  permeating  the  implanted  volume,  together  with  isolated  dislocation 
loops  [7].  Residual  damage  in  (lOO)  -oriented  substrates,  such  as  those 
used  in  Chapter  5,  is  considerably  lower  [6]. 

When  silicon  samples  are  implanted  through  an  oxide  layer,  the 
results  of  annealing  differ  significantly  from  that  of  the  bare  surface  implant. 
Well-annealed  high-dose  through-oxide  arsenic  implanted  samples  are  charac- 
terized by  stable  high  defect  density  networks  [4],  and  diffraction  patterns 
reveal  the  through-oxide  implanted  surface  layer  to  be  essentially  poly- 
crystalline in  nature. 
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An  interesting  result  from  the  earlier  studies  of  this  problem 
[4,6,7]  was  the  fact  that  implantation  of  oxygen  and  arsenic  into  bare 
silicon  surfaces  yields  exactly  the  same  type  of  disorder  seen  for  the 
through-oxide  arsenic  implanted  silicon.  The  role  of  oxygen  in  the  recovery 
of  radiation-damaged  silicon  has  been  extensively  studied.  Oxygen  is  known 
to  interact  strongly  with  irradiation-induced  defects  in  silicon  such  as  the 
lattice  vacancy  (forming  the  A-center  [77]),  and  the  As-vacancy  complex 
(the  E center)  [8].  In  fact,  it  is  the  presence  of  an  oxygen  related 
defect  in  the  recombination  spectrum  from  irradiated  silicon  that  allowed 
the  identification  of  the  oxygen  recoil  in  through-oxide  arsenic  implanted 
substrates  (Chapter  2).  With  oxygen  impurities  already  established  as  an 
important  factor  in  the  recovery  of  irradiated  silicon,  it  became  important 
to  determine  whether  or  not  the  presence  of  arsenic  was  crucial  to  the 
formation  of  the  residual  defect  structure.  If  arsenic  were  involved  in 
defect  complexes  responsible  for  the  formation  of  the  residual  defects, 
it  would  not  show  the  electrical  activation  seen  from  silicon  free  from  the 
residual  defects.  However,  the  measurements  of  the  electrically-active 
impurity  profiles  following  1000°C  annealing  (Figs.  5.6  and  5.8)  show 
arsenic  activation  for  the  through-oxide  implanted  sample  equal  to  that  seen 
from  the  bare  surface  implanted  samples.  We  must  therefore  conclude  it  is 
the  presence  of  the  excess  (recoil-implanted)  oxygen  that  is  responsible  for 
residual  defect  formation. 

To  understand  the  role  of  oxygen  in  the  formation  of  residual  defects, 
it  is  necessary  to  examine  the  behavior  of  oxygen  In  undamaged  silicon. 

Oxygen  Is  an  Impurity  commonly  found  in  Czochralski -grown  Si  [47] . The 


18  “3 

solubility  of  oxygen  in  silicon  is  approximately  2 x 10  cm  at  the  growth 
temperature  [78J.  This  solubility  decreases  with  decreasing  temperature 
until  it  reaches  a value  of  approximately  2 x lO^^cm  ^ at  1000°C  f78]. 

During  subsequent  heatings,  the  oxygen  in  solution  may  be  larger  than  the 
equilibrium  value  at  that  temperature;  if  sufficient  thermal  energy  is 
present,  the  excess  oxygen  will  then  precipitate  out  of  solution.  This 
precipitation  can  occur  in  the  form  of  Si02  precipitates  [79J,  which  would 
cause  dislocations  in  the  silicon  substrate  due  to  differences  in  physical 
properties  of  the  precipitates  from  those  of  the  surrounding  silicon 
matrix  [78].  Thus  a possible  cause  of  the  residual  damage  would  be  the 
precipitation  of  the  excess  oxygen.  Furthermiire , the  presence  of  Si02 
precipitates  would  interrupt  the  orderly  regrowth  of  the  crystalline  silicon 
from  an  amorphous  layer,  and  could  even  serve  as  nucleation  sites  for  the 
microcrys ta Is . This  hypothesis  is  clearly  consistent  with  the  fact  that 
both  transmission  electron  microscopy  and  diffraction  patterns  indicate  that 
the  recrystallized  through'oxide  implanted  layer  is  essentially  polycrysta 1 line , 
consisting  of  a collection  of  microcrystals  randomly  oriented  [4,7]. 

The  essentially  polycrystalline  nature  of  the  through-oxide  implanted 
layer  suggested  that  the  presence  of  the  excess  oxygen  affects  the  regrowth 
rate  of  the  amorphous  layer.  This  aspect  of  the  problem  was  recently  exa.mined 
by  Kennedy  e^  using  the  technique  of  helium  ion  backscattering  [80] . 

They  found  that  the  regrowth  rate  of  amorphous  silicon  decreases  exponentially 
with  increasing  oxygen  concentration,  thus  suggesting  a slight  modification 
in  the  interpretation  of  the  recovery  of  silicon  driven  amorphous  by  ion 
implantation.  Instead  of  a simple  amorphous-to-crysta lline  phase  transition, 
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epitaxial  regrowth  of  an  implanted  layer  occurs  in  competition  with  homo- 
geneous nucleation  of  the  implanted  layer  (i.e.  the  formation  of  randomly- 
oriented  silicon  crystals  at  random  locations  in  the  amorphous  layer).  For 
bare  surface  implanted  samples,  epitaxial  regrowth  apparently  predominates. 
Significantly,  the  substrate  orientation  with  the  lowest  regrowth  rate,  i.e. 
the  (ill)  direction  [76],  also  shows  the  greater  residual  disorder  after 
annealing  of  bare  surface  implanted  samples  [5] . For  samples  implanted 
through  thin  oxide  layers,  the  epitaxial  regrowth  is  slowed  by  the  presence 
of  the  recoil-implanted  oxygen,  and  thus  homogeneous  nucleation  predominates. 
Significantly,  the  residual  damage  seen  following  implantation  through  oxide 
layers  is  greatest  for  (ill)  -oriented  substrates  which  also  have  the  slowest 
regrowth  rates  for  bare  surface  implants.  Unfortunately,  it  is  not  clear 
from  these  measurements  whether  it  is  the  presence  of  SiO^  precipitates  that 
retards  the  epitaxial  regrowth  or  vd\ether  the  excess  oxygen  slows  the 
epitaxial  regrowth  rate  for  other,  unknown  reasons.  Nevertheless,  the  fact 
that  recoil  Implanted  oxygen  does  indeed  reduce  the  epitaxial  growth  rate  is 
unambigous . 

The  hypothesis  that  the  residual  damage  is  caused  by  the  slowing 
of  the  epitaxial  regrowth  rate  by  recoil  implanted  oxygen  leads  to  predictions 
that  are  verified  by  the  electrical  measurements  of  Chapter  5.  If  arsenic  is 
not  involved  in  the  defect  formation,  it  should  show  an  electrical  activation 
in  the  through-oxide  implanted  layers  equivalent  to  that  seen  for  bare  surface 
implants,  since  the  activation  of  arsenic  does  not  depend  on  the  physical 
dimensions  of  the  crystal  but  rather  on  the  concentration  of  arsenic  inside 
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a unit  volume.  Additionally,  grain  boundary  scattering  of  carriers  is 
unobservable  in  silicon  except  for  very  pure  crystals  at  very  low 
temperatures  [66,81].  For  samples  doped  to  the  high  concentrations  used 
in  this  study,  scattering  by  ionized  impurities  dominates  the  carrier 
transport  properties  of  the  crystal.  Thus,  if  the  extent  of  the  residual 
damage  were  a large  number  of  grain  boundaries,  the  carrier  mobilities 
would  be  determined  by  the  carrier  concentration  and  thus  would  yield 
mobility  profiles  identical  to  those  expected  from  bulk  silicon.  The 
experimentally-determined  arsenic  activation  and  the  mobility  profiles  are 
exactly  what  one  would  expect  by  assuming  that  the  residual  damage  is 
caused  by  random  nucleation  during  the  recrystallization  of  the  amorphous 
layer.  When  combined  with  the  channeling  measurements  of  the  epitaxial 
regrowth  rate  of  the  through-oxide  implanted  samples,  the  results  of  this 
study  lead  to  a clear  and  consistent  interpretation  of  the  cause  of  residual 
damage  in  through-oxide  implanted  silicon.  Additionally,  this  study  provides 
the  first  complete  measurement  of  the  effect  of  residual  damage  on  the 
electrical  properties  of  through-oxide  implanted  layers.  From  these  measure- 
ments it  is  seen  that  the  vast  crysta llographic  damage  created  upon  the 
annealing  of  through-oxide  implanted  silicon  has  no  measurable  effect  on 
either  the  activation  of  the  implanted  arsenic  or  the  carrier  mobilities  in 
the  implanted  layer. 

Thus,  while  this  study  has  helped  to  resolve  some  of  the  basic 
difficulties  in  understanding  the  residual  damage  caused  by  through-oxide 
arsenic  implantations,  other  questions  remain  unanswered.  For  example,  while 
grain  boundaries  and  dislocations  have  little  effect  on  carrier  mobility. 
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they  do  have  a strong  effect  on  minority  carrier  lifetimes  [82,83]. 

When  these  dislocations  occur  at  p-n  junction  boundaries,  they  serve 
as  carrier  recombination  centers  under  forward  bias  and  carrier  generation 

i 

centers  under  reverse  bias  [3],  thus  increasing  reverse  leakage.  Additionally,  i 

it  is  known  that  metallic  impurities  tend  to  precipitate  at  dislocations 
[84,85].  Such  precipitates  would  also  increase  reverse  leakage  for  p-n 
junctions.  Increased  junction  leakage  would  decrease  solar  cell  efficiency, 
reduce  the  performance  of  charge  coupled  devices,  and  lower  the  current  gain 
of  bipolar  transistors.  Since  photoluminescence  depends  strongly  on  minority 
carrier  lifetime,  the  lack  of  photoluminescence  for  anneals  above  600°C 
tends  to  indicate  that  minority  carrier  lifetimes  have  not  completely  recovered 
even  following  1000°C  annealing;  however,  these  results  should  be  supplemented 
by  other  lifetime  measurements. 

Because  dislocations  can  have  a very  strong  effect  on  minority 
carrier  lifetimes,  it  would  be  of  extreme  interest  to  examine  by  transmission 
electron  microscopy  whether  the  crystallographic  defects  caused  by  annealing 
high-dose  through-oxide  implanted  samples  propagate  into  the  substrate  following 
redistribution  by  thermal  oxidation  or  whether  this  defect  layer  is  consumed 
during  the  oxidation. 

Finally,  the  above  interpretation  is  dependent  on  the  amorphization 
of  the  surface  layer.  Thus  it  would  be  useful  to  examine  junction 
leakage  properties,  minority  carrier  lifetimes,  and  crystallographic  features 
for  through-oxide  implanted  samples  when  the  total  implanted  dose  does  not 


drive  the  substrate  amorphous. 
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For  these  reasons,  the  following  experiments  should  be  performed; 

1.  A study  should  be  made  of  the  reverse  leakage  for  arsenic -implanted 
diodes.  Results  for  diodes  defined  by  near ly-vertica 1 openings  in  the 

1 

thick  SiO^  mask  should  be  compared  with  the  results  from  diodes  defined 
by  openings  with  tapering  edges  which  would  allow  recoil  oxygen  implanta- 
tion at  the  junction  perimeter.  | 

2.  Measurements  should  be  made  of  the  minority  carrier  lifetime  for  the  | 
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5 X 10  cm  As  implants  both  into  bare  surface  and  through  Si02  to 
determine  the  effect  of  residual  damage  on  the  carrier  lifetimes  in 
these  layers . 

3.  Should  either  of  the  above  experiments  indicate  a significant  difference 
in  the  device  properties  of  the  recoil-implanted  layer,  that  experiment 
should  be  repeated  for  identical  through-oxide  implanted  samples  which 
had  been  j.irst  subjected  to  careful  gettering  to  remove  metallic 
impurities  before  implantation.  This  last  experiment  would  indicate 
whether  the  residual  damage  in  itself  was  responsible  for  the  decline 
in  device  properties,  or  whether  the  residual  damage  was  gettering  metallic 
impurities  responsible  for  the  decline  in  device  properties. 

4.  Examinations  of  the  bare  surface  and  through -oxide  implanted  samples  which 
had  been  subjected  to  redistribution  by  thermal  oxidation  should  be 
performed  by  electron  diffraction  and  by  transmission  electron  microscopy 
to  determine  whether  or  not  the  residual  damage  propagates  during  oxide 
growth . 
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5.  Examinations  should  also  be  performed  by  transmission  electron  micros- 
copy and  by  electron  diffraction  for  we  1 1 -annca led  samples  that  had 
been  implanted  tlirough  tliin  oxide  layers  to  fluences  above  and  below 
the  amorphous  limit  to  determine  if  the  sample  must  be  driven  amorphous 
for  the  residual  damage  to  occur. 

In  conclusion,  the  present  work  has  addressed  itself  to  the  origin 
and  effects  of  residual  damage  in  silicon  seen  following  through -oxide 
arsenic  implantation.  After  firmly  establishing  the  oxygen  dependence  of 
one  of  the  luminescence  peaks  seen  from  implanted  silicon,  we  have  used 
this  luminescence  to  detect  the  presence  of  recoil-implanted  oxygen  in  the 
silicon  substrate.  Electrical  measurements  were  performed  which  showed 
electrical  activation  and  carrier  .i.obility  identical  to  that  seen  from  bare 
surface  implanted  silicon.  These  measurements  were  then  found  to  be  consistent 
with  the  interpretation  that  the  residual  damage  is  caused  by  reduction  in 
the  growth  rate  of  amorphous  silicon,  favoring  the  process  of  homogeneous 
nucleation  due  to  the  presence  of  the  excess  reccil  implanted  oxygen. 
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To  obtain  the  formulas  by  which  we  convert  sheet  resistivity  and 
sheet  Hall  coefficient  data  into  carrier  concentration  and  mobility,  wo 
begin  by  considering  a structure  composed  of  two  distinct  layers,  each 
of  which  is  assumed  to  have  uniform  carrier  concentration  and  mobility 
(Figure  Al).  The  surface  layer  of  thickness  dj^  has  conductivity  and 
Hall  coefficient  Rj^2 ' application  of  a voltage  across  the  length  L 

of  this  structure  will  produce  currents  and  in  the  respective  layers, 
which  will  depend  on  the  conductivities  and  thicknesses  of  each  layer. 

If  a magnetic  field  B is  now  applied  perpendicular  to  the  sample,  each 
layer  will  develop  its  own  Hall  voltage,  which  combine  to  produce  a Hall 
voltage  for  the  entire  sample.  This  resultant  Hall  voltage  can  be  deter- 
mined by  examining  the  equivalent  circuit  of  Figure  Al.  Here  G represents 
the  conductance  of  each  layer  in  the  direction  parallel  to  the  contacts 


Ld  Ld 

S = °1  — S = ^^2  — 


(Al) 


From  Kirchoff's  laws,  the  equivalent  voltage  V seen  at  the  contacts  is  then 
given  by 


V = 


Gi+G^ 


(A2) 


where  the  Hall  voltages  are  given  by 


V bW 

^ ■ '^l 


2 d„ 


(A3) 
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Figure  A.l.  Top:  Hypothetical  two-layer  structure  used  to  determine 

sheet  resistivity  and  sheet  Hail  coefficient  for  composite 
sample . 

Bottom;  Equivalent  circuit  used  to  determine  Hall  voltage 
for  two-layer  structure  shown  above. 


A 


r 


(^1^1  ■'•^2^2^' 


(A  10) 


Since  layers  one  and  two  are  in  parallel  in  the  direction  of  current  flow 
their  conductivities  add  to  give  the  conductivity  of  the  entire  sample 

adL  c^d^L 

V = (All) 


CT  = - (a^d^i-a^d^) 


(A12) 


with  the  conductance  of  the  entire  sample. 

Generalizing  the  results  (AlO)  and  (A12)  to  a total  of  j independent 
layers,  we  obtain 


X=1  1 


‘ j 2 


(A  13) 


ff  = 7 S a.d.  ( 

d . , 1 1 ' 

1=1 

If  we  now  assume  that  the  Hall  mobility  u is  equal  to  the  conductivity 

H 

mobility  a and  also  assume  that  the  currents  is  each  layer  are  carried  by 
majority  carriers  only,  we  then  have 


<J.  = n.qu.  and  R = 

111  H,  n.q 


(A  14) 


(A  15) 
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so  that 


a = — T.  n.qn..d.  and  R„  = d 
d . , 1 1 1 H 

1=1 


r n.qa.d. 
i=l  ^ ^ ^ 


( S n.qu.d. ) 

Li  = i ^ ^ 


- 16 

with  q equal  to  the  electronic  change  (1.602x10  coul).  The  sheet 


(A16) 


resistivity  p for  our  structure  of  j layers  is  defined  by 


-1 


J 


P = (ad)  = ( ^ n qu  ) 

S . .,11 

J 1=1 


-1 


(A  17) 


and  the  sheet  Hall  coefficient  for  our  structure  of  j layers  is 


R 


J 2 

S n.qu.d. 

_ ^ ^ 


s d 3 , 

^ ( r n.qu.d. ) 

. 1 1 1 
1=1 


(A18) 


By  examining  the  forms  of  the  above  expressions  we  can  see  that 


and 


r J 1 - 

S n.qu.d. 

i=l  1 1 ^ 

1 

1 s2 

J -1 

( r n.qu.d. ) 

( S n. qu. d . ) 

Li=i  ^ ^ _ 

Li=i  111  J 

(A19) 


1 


D 


S . 

J 


j 

S n. qU. d . 
i=l  " " " 


(A20) 


From  the  last  two  expressions  we  see  that  if  we  measure  both  the  sheet 
resistivity  and  sheet  Hall  coefficient  for  our  structure  consisting  of  j 
hypothetically  distinct  layers  and  then  remove  one  layer,  our  structure  will 
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now  consist  of  j-1  layers  and  have  sheet  resistivities  and  sheet  Hal 
coefficients  given  by 


5 = 7]  n .qu.d. 

s . 1 1 1 

J-1  1=1 


(A21) 


^ 2 
I!  n.qu.d, 

i=l  " ^ " 


j-1  j-1 


( y.  n.qu.d.  ) 
. ,111 
1=1 


(A22) 


so  that  the  changes  in  the  quantities  defined  in  (A19)  and  (A20)  will  become 


A(R  /o  ). 
s s j 


R R 
s j 

2 


■j-1 


= y n.qu.d.  - y n.qu.d.  = n.qu.d.  (A23) 


s . 
J 


^j-1 


i=l 


1 11 


i=l 


J J J 


and 


A(l/o  ).  = — 

s 2 0.  0 


j j-1 

= ? n.qu.d.  - y n.qu.d.  = n.qu.d.  (A24) 

. 1 1 1 1 .,111  J J J 

s.  s.  . 1=1  1=1  J J J 

J J-1 


Thus,  we  obi.ain  the  desired  results  by  considering  dividing  (A23)  by  (A24) 
and  then  using  this  result  in  (A24) 


and 


u. 

J 


A(r  /o^) 
s s 


j 


A(1/d). 


A(l/0g) 

qdu . 


(A25) 


n . 
J 


(A26) 
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